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Abstract 

Future  missions  to  Mars  may  require  pin-point  landing  precision,  possibly  on  the 
order  of  tens  of  meters.  The  ability  to  reach  a  target  while  meeting  a  dynamic 
pressure  constraint  to  ensure  safe  parachute  deployment  is  complicated  at  Mars  by 
low  atmospheric  density,  high  atmospheric  uncertainty,  and  the  desire  to  employ  only 
bank  angle  control.  The  vehicle  aerodynamic  performance  requirements  and  guidance 
necessary  for  a  0.5  to  1.5  lift-to-drag  ratio  vehicle  to  maximize  the  achievable  footprint 
while  meeting  the  constraints  are  examined. 

A  parametric  study  of  the  various  factors  related  to  entry  vehicle  performance  in 
the  Mars  environment  is  undertaken  to  develop  general  vehicle  aerodynamic  design 
requirements.  The  combination  of  low  lift-to-drag  ratio  and  low  atmospheric  density 
at  Mars  results  in  a  large  phugoid  motion  involving  the  dynamic  pressure  which 
complicates  trajectory  control.  Vehicle  ballistic  coefficient  is  demonstrated  to  be  the 
predominant  characteristic  affecting  final  dynamic  pressure.  Additionally,  a  speed 
brake  is  shown  to  be  ineffective  at  reducing  the  final  dynamic  pressure. 

An  adaptive  precision  entry  atmospheric  guidance  scheme  is  presented.  The  guid¬ 
ance  uses  a  numeric  predictor-corrector  algorithm  to  control  downrange,  an  azimuth 
controller  to  govern  crossrange,  and  an  analytic  control  law  to  reduce  the  final  dy¬ 
namic  pressure.  Guidance  performance  is  tested  against  a  variety  of  dispersions,  and 
the  results  from  selected  test  cases  are  presented.  Precision  entry  using  bank  angle 
control  only  is  demonstrated  to  be  feasible  at  Mars.  '' 
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Chapter  1 
Introduction 


Future  exploratory  missions  to  the  surface  of  Mars  will  require  precise  control  of  the 
terminal  landing  conditions.  Current  studies  for  a  robotic  mission,  such  as  the  Mars 
Rover  Sample  Return  (MRSR)  study  [1],  levy  requirements  of  100  meter  accuracy 
in  order  to  land  in  geologically  interesting  areas  or  to  avoid  hazards  that  have  been 
detected  from  orbit.  Later  missions  to  Mars  will  most  likely  involve  the  emplacement 
of  manned  habitats  and  their  subsequent  resupply.  The  landing  accuracy  required 
for  these  missions  will  be  on  the  order  of  30  meters  [2]. 

Precision  entry  and  landing,  such  as  the  Space  Shuttle  frequently  performs  at 
Earth,  has  never  been  done  at  Mars.  Relatively  little  work  has  been  done  to  define 
vehicle  and  guidance  requirements  for  such  a  mission.  The  task  is  complicated  at  Mars 
by  a  less  dense  atmosphere,  the  greater  uncertainty  in  environmental  conditions,  and 
the  need  for  a  high  degree  of  autonomy.  This  thesis  attempts  to  define  basic  vehicle 
performance  requirements  for  achieving  entry  and  precision  landing,  and  then  presents 
a  guidance  algorithm  that  maximizes  the  achievable  footprint  for  a  given  vehicle  while 
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minimizing  the  landing  errors. 

In  support  of  these  goals,  this  chapter  presents  background  regarding  entry  and 
formulates  the  entry  problem.  Section  1.1  discusses  the  different  entry  options  avail¬ 
able  and  provides  an  overview  of  the  general  entry  mission.  Given  a  specific  mission, 
Section  1.2  then  introduces  the  constraints  and  requirements  which  define  a  successful 
entry.  There  are  several  possible  methods  of  controlling  the  vehicle  in  order  to  achieve 
these  mission  requirements.  Section  1.3  discusses  these  possibilities  and  presents  the 
chosen  method  of  control.  With  the  specific  entry  constraints  and  controls  defined,  the 
approach  used  to  solve  the  entry  problem  and  an  overview  of  the  thesis  are  contained 
in  Section  1.4. 

1.1  Mission  Overview 

The  exploration  of  Mars  has  been  of  interest  to  man  for  many  years.  Current  interest 
revolves  around  the  human  exploration  of  Mars,  and  President  Bush’s  Lunar/ Mars 
Exploration  Initiative  provides  the  framework  by  which  this  will  take  place.  The 
plan  calls  for  a  series  of  missions  which  begins  with  unmanned  probes,  leads  to  the 
initial  manned  landing,  and  eventually  achieves  a  permanent  presence.  However, 
many  technological  issues  must  be  addressed  before  these  missions  become  reality. 
Precision  entry  and  larding  is  one  of  these  issues. 

Because  little  work  has  been  done  in  this  arena,  there  are  no  definite  mission  pro¬ 
files  to  work  with.  In  general,  there  are  two  ways  that  a  vehicle  on  an  interplanetary 
trajectory  to  Mars  can  reach  the  surface.  The  vehicle  can  enter  directly  from  the  hy- 


22 


perbolic  approach  orbit,  or  it  can  capture  into  a  lower  energy  orbit  and  then  deorbit 
to  the  surface.  The  capture  approach  uses  more  fuel,  but  it  has  several  advantages. 
Besides  reducing  the  entry  velocity  which  lowers  the  heating  rates  and  g-loads,  this 
approach  allows  time  for  on-orbit  determination  of  a  safe  landing  site.  This  may  be 
necessary  for  the  early  missions.  Once  landing  sites  and  base  camps  are  established 
later  in  the  -'rogram,  the  capture  approach  may  still  be  desirable  based  on  mission 
planning  considerations. 

Once  the  vehicle  reaches  the  atmosphere  of  Mars,  there  are  two  general  classes 
of  entry  trajectories  which  can  be  flown.  The  trajectory  can  be  lifting  (employing  a 
vehicle  which  creates  aerodynamic  lift)  or  non-lifting  (using  a  vehicle  which  develops 
no  lift).  The  non- lifting  trajectory  is  generally  shorter,  subject  to  smaller  dispersions, 
and  thus  more  accurate.  The  original  Minuteman  Intercontinental  Ballistic  Missile 
warheads  used  this  type  of  trajectory.  However,  the  non-lifting  trajectory  results 
in  entry  and  terminal  conditions  which  are  very  severe,  and  it  is  not  feasible  for 
manned  flight.  The  lifting  entry  trajectory,  as  used  in  the  Apollo,  Shuttle,  and  Viking 
programs,  generally  moderates  the  entry  conditions  at  the  cost  of  being  more  sensitive 
to  dispersions.  A  larger  lifting  capacity  allows  greater  reductions  in  the  g-loads  and 
final  velocities  by  increasing  control  over  the  trajectory. 

Although  entry  at  Mars  has  been  successfully  demonstrated  by  the  American 
Viking  landers  [3]  and  the  Soviet  probe  Mars  6  [4j,  these  entries  were  not  precision 
entries  by  any  means.  The  3cr  landing  error  ellipse  for  the  low-lift,  open-loop  guidance, 
Viking  1  trajectory  was  on  the  order  of  120  km  x  60  km  [3].  In  order  to  provide  a  low 


23 


terminal  velocity,  allow  for  precise  targeting,  and  provide  a  more  benign  environment 
during  entry,  a  guided  mid-to-high  lift-to-drag  ratio  (L/D)  entry  vehicle  with  the 
capability  to  significantly  shape  the  trajectory  is  desired. 


1.2  Entry  Requirements 

One  landing  scenario  studied  under  the  MRSR  [1]  mission  includes  a  three  phase 
profile  for  reducing  the  vehicle  velocity  while  actively  guiding  to  a  landing  site.  This 
is  the  scenario  assumed  in  this  study,  and  the  three  successive  phases  are, 

•  Entry  Phase 

•  Parachute  Phase 

•  Powered  Descent  Phase 

The  entry  phase  begins  once  atmospheric  forces  start  to  affect  the  trajectory.  Dur¬ 
ing  this  stage,  the  vehicle  uses  its  aerodynamic  capabilities  to  fly  towards  the  target 
while  decelerating  from  orbital  to  near  subsonic  velocities.  The  entry  phase  ends 
when  the  vehicle  deploys  its  parachute  at  some  predetermined  altitude  or  combina¬ 
tion  of  altitude  and  velocity.  The  parachute  phase  is  very  short  compared  to  the 
approximately  one-half  hour  entry  phase  and  lasts  only  about  a  minute.  The  purpose 
of  the  parachute  phase  is  to  further  reduce  the  vehicle  velocity.  After  the  parachute  is 
jettisoned,  the  descent  engines  are  turned  on,  and  the  powered  descent  phase  begins. 
During  this  phase,  the  descent  engines  remove  the  vertical  velocity  from  the  trajec¬ 
tory  and  maneuver  the  vehicle  horizontally  to  the  landing  site.  This  study  focuses 


24 


on  vehicle  capabilities  and  guidance  performance  during  the  first  phase — the  entry 
stage. 

A  set  of  entry  conditions  from  [5]  is  used  to  define  the  requirements  in  this  study. 
These  requirements  are  summarized  in  Table  1.1  and  explained  in  the  following  para¬ 
graphs.  These  requirements,  composed  of  final  conditions  and  continuous  constraints, 
are  subject  to  change  and,  at  this  point  in  mission  development,  should  be  regarded 
only  as  typical  guidelines. 


Table  1.1.  General  Performance  Requirements 


Terminal  Conditions 
Down  Range  Error  <  1  km 
Cross  Range  Error  <  1  km 
Final  Dynamic  Pressure  <  1000  Pascals 
Final  Flight  Path  Angle  between  —20°  and  0° 

Continuous  Constraints 
G- loads  <  3  Earth  g’s 
Low  Heating  Rate 
Large  Crossrange  Capability 


The  landing  accuracy  expected  of  the  complete  entry  and  landing  profile  varies, 
depending  upon  the  specific  mission  requirements,  and  may  be  as  small  as  30  m  [2]. 
This  pin-point  landing  capability  drives  guidance  accuracy  during  the  entry  phase. 
Since  dispersions  due  to  wind  are  expected  during  the  parachute  phase,  and  current 
plans  call  for  less  than  2  km  of  lateral  maneuvering  capability  during  the  powered 
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descent  phase  [1],  the  entry  phase  must  guide  the  vehicle  to  within  approximately  1 
km  of  the  desired  parachute  deployment  location. 

The  final  dynamic  pressure  limitation  ensures  the  structural  integrity  of  the  para¬ 
chute  and  limits  the  opening  shock  experienced  by  the  vehicle.  The  parachute  is  also 
a  driver  behind  the  flight  path  angle  requirement.  The  entry  guidance  must  leave 
the  vehicle  on  a  nearly  horizontal  trajectory  in  order  to  minimize  variations  in  actual 
parachute  deployment  altitude  and  allow  sufficient  time  for  parachute  inflation  and 
vehicle  deceleration. 

The  g-load  limitation  is  intended  to  provide  comfortable  entry  conditions  for  as¬ 
tronauts  who  might  be  in  a  weakened  state  following  a  long  journey  from  Earth,  as 
well  as  to  maintain  vehicle  integrity.  The  heating  rate  must  be  kept  low  in  order  to 
reduce  thermal  protection  subsystem  requirements. 

An  additional  requirement  for  entry  is  a  large  crossrange  capability.  A  crossrange 
capability  is  necessary  to  reach  landing  sites  at  a  higher  latitude  than  the  inclination 
of  the  orbit,  or  to  reach  landing  sites  not  in  the  orbital  plane  at  the  time  of  deorbit. 
A  greater  crossrange  capability  increases  the  number  of  landing  sites  which  can  be 
reached  on  a  given  orbital  pass.  This  increases  the  number  of  deorbit  opportunities 
to  a  given  landing  site  and  enhances  mission  flexibility.  Downrange  capability  is  not 
included  as  a  constraint  because  all  vehicles  have  some  inherent  downrange  ability 
and  downrange  targeting  is  easily  changed  by  delaying  the  time  of  the  deorbit  burn. 
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1.3  Trajectory  Control 

The  two  major  forces  which  act  on  the  vehicle  during  entry  are  gravitational  and 
aerodynamic.  In  order  to  modify  the  trajectory  and  meet  the  constraints,  the  entry 
vehicle  must  be  able  to  modulate  the  two  components  of  the  aerodynamic  force,  lift 
and  drag. 

Figure  1.1  from  [6]  displays  a  general  range  of  possible  entry  vehicles  and  their 
relative  lift-to-drag  ratio,  L/D,  and  ballistic  coefficient,  Cb ■  The  higher  L/D  vehicles 
are  desirable  because  they  provide  a  higher  degree  of  trajectory  control,  but  they 
have  sharp  leading  edges  or  wings  for  which  there  are  substantial  component  packag¬ 
ing  and  mass  penalties.  The  lower  L/D  vehicles  have  excellent  mass  and  packaging 
characteristics,  but  their  spherical  and  simple  cone  shapes  do  not  generate  enough  lift 
to  greatly  modify  the  trajectory.  In  between  these  vehicles  is  the  general  class  which 
is  being  considered  for  Mars  entry — the  lifting  biconics. 

A  biconic  vehicle,  as  shown  in  Figure  1.2  [7],  is  axially  symmetric  and  generates 
lift  by  the  displacement  of  the  center  of  gravity  from  the  axis  of  symmetry.  The 
amount  of  lift  generated  can  be  controlled  by  varying  the  angle  of  attack,  a,  and 
this  can  be  accomplished  through  the  use  of  a  body  flap,  reaction  control  jets  or  by 
movement  of  the  center  of  gravity.  The  body  flap  adds  weight  and  complexity,  while 
the  reaction  control  jet  solution  is  expensive  from  a  fuel  standpoint.  Additionally, 
the  shifting  of  the  center  of  gravity  is  difficult  to  control.  This  makes  it  desirable  to 
avoid  controlling  the  magnitude  of  the  lift  with  the  angle  of  attack  and  to  simply  fly 
at  the  constant  trim  angle  of  attack,  at,  determined  by  the  placement  of  the  center 
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HYPERSONIC  LIFT-TO-DRAG  RATIO  (UD) 


HYPERSONIC  BALLISTIC  COEFFICIENT  IM/CpA) 


Figure  1.1.  Range  of  Entry  Vehicle  Configurations 
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Figure  1.2.  Typical  Biconic  Vehicle 

of  gravity. 

If  it  is  decided  not  to  directly  control  the  amount  of  lift  being  generated  by  mod¬ 
ulating  a,  it  is  still  possible  to  control  the  trajectory  of  a  biconic  vehicle  by  simply 
varying  the  direction  of  the  lift  created.  This  is  done  by  rotating  the  vehicle,  and  thus 
the  lift  vector,  about  the  atmosphere-relative  velocity  vector,  using  reaction  control 
jets  or  aerodynamic  flaps.  Reaction  jets  are  the  preferred  method  of  controlling  bank 
angle  because  they  are  effective  when  the  aerodynamic  forces  are  too  low  to  gener¬ 
ate  the  necessary  rates,  and  because  they  will  probably  already  be  on  the  vehicle 
to  provide  control  on-orbit.  Thus,  using  only  reaction  control  jets  avoids  an  extra 
system. 

The  second  component  of  the  aerodynamic  force,  drag,  can  be  directly  affected 
by  modulating  the  vehicle  surface  area  or  the  coefficient  of  drag.  A  speed  brake 
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can  modify  both  of  these,  but  at  the  cost  of  increased  complexity.  Drag  control  can 
be  indirectly  achieved  by  using  the  bank  angle  to  vary  the  vertical  lift,  and  thus 
the  altitude  and  density  at  which  the  vehicle  is  flying.  This  form  of  drag  control  is 
desirable  because  it  does  not  require  structural  additions  to  the  vehicle. 

Bank  angle  control  alone  allows  control  authority  over  both  lift  and  drag  during 
entry,  without  adding  undue  weight  or  complexity  to  the  vehicle  by  requiring  angle  of 
attack  control.  Therefore,  for  the  purposes  of  this  report,  the  entry  vehicle  is  assumed 
to  be  a  biconic  vehicle  with  bank  angle  control  only.  The  extent  of  this  authority  and 
whether  it  is  sufficient  for  precision  landing  is  to  be  analyzed  in  this  study. 

1.4  Thesis  Overview 

In  analyzing  the  precision  entry  and  landing  problem,  this  thesis  uses  simulations  of 
entry  at  Mars.  The  simulator  and  the  models  used  to  perform  the  presented  analyses 
are  described  in  Chapter  2. 

After  the  development  of  the  simulator,  the  thesis  investigates  two  aspects  of  the 
precision  entry  and  landing  problem.  Tfe  first  aspect  studied  is  the  performance  of 
entry  vehicles  in  the  Martian  atmosphere.  The  ranging  and  maneuvering  capabilities 
and  final  conditions,  which  are  a  result  of  various  vehicle  configurations,  are  analyzed 
in  this  portion  of  the  thesis.  From  the  performance,  basic  vehicle  design  require¬ 
ments  will  be  defined.  Chapter  3  examines  entry  vehicle  performance  at  Mars  by 
parametrically  varying  vehicle  and  environment  factors.  The  important  parameters 
for  achieving  specific  performance  objectives  are  identified  here,  and  this  part  serves 
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as  an  independent  collection  of  knowledge  concerning  entry  at  Mars. 

The  second  aspect  of  this  thesis  concerns  the  ability  of  a  given  entry  vehicle  to 
reach  a  target  at  Mars  without  violating  the  constraints.  This  part  consists  of  the 
development  and  testing  of  a  guidance  algorithm  for  achieving  precision  entry  at  Mars. 
Chapter  4  describes  previous  research  relevant  to  the  guidance  and  then  presents  its 
design.  The  algorithm  attempts  to  maximize  the  landing  footprint  for  a  given  vehicle 
while  controlling  downrange,  crossrange  and  final  dynamic  pressure.  In  Chapter  5, 
the  performance  of  the  guidance  is  stress-tested  against  severe  dispersions  in  order 
to  demonstrate  the  capabilities  and  disadvantages  of  using  only  bank  angle  control 
during  entry. 

The  conclusions  drawn  from  both  parts  of  this  thesis  are  summarized  in  Chapter  6. 
The  lessons  learned  from  this  study  are  presented,  and  topics  for  future  research  are 
examined. 
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Chapter  2 


Simulation 

2.1  Introduction 

This  chapter  describes  the  computer  simulation  and  associated  models  used  to  per¬ 
form  the  analyses  presented  in  this  thesis.  Open-loop  and  closed-loop  guided  entry 
trajectories  were  generated  using  a  three  degree-of-freedom  (3-DOF)  digital  simu¬ 
lator.  The  simulator  was  developed  at  The  Charles  Stark  Draper  Laboratory,  Inc. 
using  the  MAC  programming  language  on  an  IBM  3090  mainframe  computer.  Using 
a  vehicle  model  and  an  environment  model,  the  equations  of  motion  are  integrated  in 
time  from  the  prescribed  initial  conditions  down  to  a  designated  terminal  condition. 
The  simulator  uses  a  4th-order  Runge-Kutta  integrator  with  a  one  second  time  step. 

The  simulator  results  are  dependent  upon  the  specific  models  used,  and  these 
models  are  described  in  the  remainder  of  the  chapter.  Section  2.2  presents  the  Martian 
environment  models  used.  Included  in  Section  2.2  are  the  physical  characteristics  of 
the  planet,  its  atmosphere,  and  the  expected  dispersions.  Descriptions  of  the  two 
different  methods  of  vehicle  definition  used  and  the  associated  vehicle  dispersions  are 
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contained  in  Section  2.3. 


2.2  Mars  Environment  Models 

The  important  parameters  modeled  in  the  environment  are  planet  shape,  gravity, 
atmospheric  density,  and  winds.  In  the  simulator,  the  vehicle  altitude  above  the 
surface  is  determined  relative  to  a  reference  ellipsoid  which  is  defined  by  equatorial 
and  polar  planet  radii.  This  same  reference  ellipsoid  is  used  to  determine  a  target 
vector  based  upon  the  landing  site  location  on  the  planet  surface.  However,  for 
the  purpose  of  defining  the  final  surface  range  between  the  vehicle  and  target,  a 
spherical  planet  model  is  used.  This  is  done  for  the  sake  of  simplicity,  and  the  error 
introduced  is  small  because  the  planet  is  approximately  circular  for  the  short  ranges 
considered.  The  expression  which  is  used  to  model  gravitational  acceleration  reflects 
the  non-spherical  nature  of  Mars  by  including  J2  effects.  The  value  of  the  gravitational 
constant,  /xm,  and  other  important  physical  constants  used  in  the  simulator  are  found 
in  Table  2.1. 

The  atmospheric  density  at  Mars  varies  with  season,  latitude  and  time  of  day, 
in  addition  to  altitude.  Despite  these  variations,  it  is  important  to  have  a  standard 
density  atmosphere  with  respect  to  altitude  to  serve  as  a  reference.  The  simulator 
uses  the  revised  Committee  on  Space  Research  (COSPAR)  Northern  Hemisphere 
mean  Mars  atmosphere,  described  in  Reference  [8],  as  a  nominal  density  model.  It  is 
implemented  as  a  table  look-up.  The  Martian  atmosphere  is  significantly  less  dense 
than  Earth’s  atmosphere.  As  can  be  seen  in  the  comparison  between  the  Maxs  mean 


34 


Table  2.1.  Mars  Environment  Constants 


Gravitational 

Parameter 

Vm  = 

4.28282804  xlO13 

TO3 

Gravitational 

Radius 

r9  = 

3393.4 

km 

Equatorial 

Radius 

req  — 

3393.4 

km 

Polar 

Radius 

r?l  = 

3375.8 

km 

Rotation 

Rate 

“>m  = 

2t 

88643 

rad 

sec 

Gravitational 

Perturbation 

h  = 

0.001965 

model  and  a  standard  Earth  model  in  Figure  2.1  [8],  Mars’  surface  density  is  only  1% 
of  Earth’s  surface  density.  This  lower  density  decreases  the  amount  of  aerodynamic 
forces  produced  by  an  entry  vehicle  at  a  given  altitude  and  velocity,  and  thus  reduces 
the  ability  of  the  vehicle  to  modify  its  trajectory  to  reach  a  target.  Additionally,  the 
decrease  in  drag  causes  difficulty  in  getting  the  vehicle  to  slow  down. 

In  addition  to  planning  for  the  lower  density  at  Mars,  entry  designers  must  also 
deal  with  the  lack  of  certainty  concerning  the  actual  density  which  will  be  encountered. 
Although  data  returned  from  Viking  1  and  2  and  Mars  6  has  greatly  reduced  the 
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Figure  2.1.  Comparison  of  Mars  and  Earth  Atmospheric  Densities 
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initial  uncertainty  in  the  density,  these  three  missions  constitute  a  limited  data  base, 
and  the  actual  density  could  significantly  vary  from  these  estimates.  Because  of 
this,  it  is  important  to  test  vehicle  and  guidance  performance  over  a  large  range  of 
dispersed  atmospheres.  Additionally,  the  atmospheric  density  can  experience  large 
and  unpredictable  short-term  variations  from  nominal  value.  It  is  important  to  test 
for  these  cases,  too.  The  different  profiles  which  will  be  used  to  analyze  performance 
in  off-nominal  atmospheres  are: 

1.  Revised  COSPAR  Cool-Low  Model 

2.  Revised  COSPAR  Warm-High  Model 

3.  Viking  1  Entry  Trajectory  Data 

4.  Viking  2  Entry  Trajectory  Data 

5.  Constant  Density  Bias — 50%  Thin  Northern  Hemisphere  Mean  Model 

6.  Constant  Density  Bias — 100%  Thick  Northern  Hemisphere  Mean  Model 

7.  Early  Density  Shear — 25%  Thin  Northern  Hemisphere  Mean  Model 

8.  Late  Density  Shear — 10%  Thin  Northern  Hemisphere  Mean  Model 

The  deviation  of  the  Cool-Low  and  Warm-High  models  from  the  nominal  density 
is  shown  in  Figure  2.2,  while  the  deviation  of  the  Viking  1  and  2  profiles  from  nominal 
is  contained  i:  Figures  2.3  and  2.4.  All  of  these  figures  are  from  Reference  [8]. 

The  density  biases  and  shears  are  based  on  the  Northern  Hemisphere  mean  model. 
The  early  shear  is  defined  by  a  25%  thin  atmosphere  from  entry  interface  until  50 
kilometers  altitude — at  which  time  the  density  abruptly  returns  to  nominal.  The  late 
shear  consists  of  a  nominal  atmosphere  until  25  kilometers  altitude  where  the  density 
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Figure  2.3.  Viking  1  Profile  Density  Variation 
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Figure  2.4.  Viking  2  Profile  Density  Variation 
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becomes  10%  thin  and  remains  so  down  to  the  surface.  The  late  shear  is  smaller 
percentage-wise  than  the  early  shear  because  larger  shears  are  less  likely  at  lower 
altitudes.  However,  the  late  shear  has  the  potential  to  cause  greater  errors  because 
the  vehicle  has  little  time  remaining  to  correct  for  this  dispersion. 

For  the  nominal  entry  case,  the  entire  atmosphere  is  assumed  to  be  rotating  uni¬ 
formly  with  Mars,  and  no  planet-relative  winds  are  included.  However,  winds  can  be 
an  important  factor  in  vehicle  performance,  and  dispersed  wind  test  cases  consisting 
of  a  constant  50  m/s  wind  from  each  of  the  four  cardinal  directions  are  used.  The  50 
m/s  value  was  chosen  as  a  reasonable  value  compromising  between  the  high  speed  jet 
stream  winds  which  the  vehicle  experiences  for  a  brief  period  of  time  and  the  lower 
speed  winds  generally  found  elsewhere.  Reference  [8]  indicates  that  for  every  season 
on  Mars,  the  surface  winds  recorded  by  the  Viking  landers  were  less  than  17  m/s 
ninety-nine  percent  of  the  time.  Error  analysis  of  the  Viking  1  trajectory  data  indi¬ 
cates  that  the  error  due  to  wind  is  attributable  to  the  equivalent  of  a  constant  30  m/s 
wind  from  the  east  [3]. 

2.3  Entry  Vehicle  Models 

Because  there  are  two  different  studies  described  in  this  report — one  parametri¬ 
cally  exploring  vehicle  performance  capabilities  and  the  other  examining  guidance 
performance — it  was  desirable  to  employ  two  separate  methods  of  defining  the  entry 
vehicle.  One  method  defines  a  generic  biconic  lifting  entry  vehicle  based  solely  on 
constant  values  of  lift-to-drag  ratio  and  ballistic  coefficient.  The  other  method  uses 
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vehicle  mass,  surface  reference  area,  and  the  aerodynamic  coefficients,  Cl  and  Cd-i 
a  function  of  angle  of  attack  for  a  specific  configuration  to  distinguish  the  vehicle. 

L/D  and  CB  fully  define  the  vehicle  aerodynamic  characteristics  for  the  trans¬ 
lational  equations  of  motion  and  are  therefore  the  most  important  parameters  in 
assessing  vehicle  capability.  The  first  method  of  vehicle  definition  is  useful  because 
it  facilitates  the  parametric  examination  of  different  vehicle  classes  without  access 
to  specific  configuration  data.  The  general  entry  trajectory  studies  in  Chapter  3 
are  performed  using  this  generic  vehicle.  The  second  method  requires  more  specific 
configuration  information,  but  is  useful  in  examining  the  performance  that  can  be 
extracted  from  a  particular  vehicle.  This  method  of  vehicle  definition  is  used  later 
in  Chapter  3  for  just  this  purpose.  The  guidance  performance  analysis  in  Chapter  5 
also  uses  this  method  of  vehicle  definition. 

In  order  to  realistically  assess  the  guidance  performance,  the  bank  angle  control 
response  was  modeled  as  a  phase  plane  having  bank  rate  and  acceleration  limits.  In 
this  model,  the  maximum  roll  rate  was  20  degrees  per  second,  and  the  maximum  roll 
acceleration  was  5  degrees  per  second  squared. 

Although  the  aerodynamic  characteristics  of  the  entry  vehicle  will  be  analytically 
and  empirically  predicted,  the  small  errors  expected  in  this  knowledge  can  cause 
significant  performance  differences.  The  Viking  lander  axial  force  coefficients  were 
only  in  error  by  approximately  1%,  and  this  resulted  in  a  6  to  10  km  range  error  [3]. 
If  the  Viking  trajectories  had  been  longer,  as  will  future  precision  landing  trajectories, 
the  range  error  would  have  grown  much  larger.  Since  the  entry  flight  regime  is  beyond 
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that  of  current  wind  tunnels,  and  computational  fluid  dynamics  research  has  not 
sufficiently  advanced  to  allow  accurate  prediction,  there  will  be  errors.  There  will  also 
be  a  large  uncertainty  in  the  trim  angle  of  attack,  and  this  will  induce  a  significant 
known  L/D  dispersion.  Therefore,  ±  10%  dispersions  in  lift  and  drag  coefficients  will 
be  considered. 

The  vehicle  heating  rate  during  entry  is  modeled  by  an  empirical  formula  found 
in  [9].  The  equation  estimates  the  stagnation  heating  rate  on  the  nose  of  the  entry 
vehicle,  in  W/sq  cm,  as, 


A  542.01 

Q~ 


P_ 

Pile 


Vrel  \3'5 

3048/ 


(2.1) 


where  Rpt  is  the  vehicle  nose  radius  in  units  of  meters  and  vTei  is  in  units  of  m/s.  The 
values  used  in  the  simulator  for  the  constants  in  the  equation  are, 

Rn  =  0.3048  m 
and 


p,ie  =  1.226  kg/m3 

The  integrated  heat  load  during  entry  is  then  computed  as  the  time  integral  of  the 
heating  rate  equation  as, 


Q{t)=  f  Q(t)<1t 

Jo 


(2.2) 


43 


44 


Chapter  3 

Entry  Vehicle  Capabilities 


3.1  Introduction 

The  ability  to  land  with  pin-point  precision  is  highly  dependent  on  the  entry  vehicle 
capabilities.  This  chapter  examines  lifting  vehicle  entry  trajectories  at  Mars  and 
looks  at  the  effect  of  different  vehicle  and  environmental  factors  on  performance. 
Section  3.2  describes  how  the  studies  in  this  chapter  were  performed.  Section  3.3 
examines  the  entry  ranging  and  final  dynamic  pressure  performance  over  a  range 
of  vehicle  configurations,  while  Section  3.4  investigates  the  same  performance  for  a 
specific  vehicle.  The  last  section  examines  the  ability  of  a  vehicle  to  divert  late  in  the 
Mars  entry  trajectory  in  order  to  reach  an  updated  target  location. 

3.2  Method  of  Performance  Analysis 

The  performance  criteria  used  in  these  evaluations  are  downrange  and  crossrange¬ 
ranging  capability-  and  final  dynamic  pressure  because  they  are  the  primary  entry 
constraints  and  will  be  used  as  a  basis  of  guidance  design  requirements.  Downrange  is 
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Table  3.1.  Entry  Interface  Conditions 


Altitude 

125.0 

km 

Inertial  Velocity 

3550.0 

m/s 

Flight  Path  Angle 

-3.3 

deg 

Latitude 

0.0 

deg 

Longitude 

0.0 

deg 

Heading 

90.0 

deg 

Vacuum  Apoapsis 

500.0 

km 

Vacuum  Periapsis 

5.0 

km 

measured  from  the  initial  position  vector  to  the  projection  of  the  final  position  vector 
on  the  initial  orbital  plane.  Crossrange  is  measured  as  the  perpendicular  distance  of 
the  final  position  vector  from  the  same  orbital  plane.  All  of  these  vectors  are  specified 
in  a  planet-fixed  coordinate  system. 

Next,  the  initial  and  final  conditions  must  be  selected.  For  the  purposes  of  this 
study,  the  vehicle  is  initially  in  a  500  km  circular  orbit  and  deorbits  by  targeting  for 
a  vacuum  periapsis  altitude  of  5  km.  This  results  in  a  flight  path  angle  of  -3.3°and  a 
velocity  of  3550  m/s  at  entry  interface  (125  km  altitude).  For  simplicity,  the  vehicle 
location  is  0°  longitude  and  0°  latitude  at  entry  interface,  and  the  vehicle  heading  is 
due  east.  These  nominal  initial  conditions,  summarized  in  Table  3.1,  are  used  for  all 
test  cases  in  this  chapter  with  the  exception  of  those  dispersed  in  initial  flight  path 
angle. 

Since  the  aeromanuevering  trajectory  is  defined  to  end  at  parachute  deployment, 
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this  is  where  the  terminal  conditions  are  measured.  For  this  study,  the  parachute 
is  deployed  at  5  km  altitude,  and  this  is  chosen  to  be  the  trajectory  termination 
condition. 

Now  that  the  initial  and  final  conditions  have  been  specified,  a  common  control 
profile  which  demonstrates  the  full  capabilities  of  the  vehicle  must  be  chosen.  Flying 
constant  full  lift-up  will  provide  the  maxirr.  im  downrange  possible  for  a  given  vehicle, 
but  trajectories  with  a  constant  bank  profile  do  not  achieve  maximum  crossrange 
performance.  A  vehicle  flying  a  constant  full  lift  out-of-plane  (0  =  90°)  bank  profile 
will  turn  quickly,  but  will  fail  to  go  far  in  crossrange  because  there  is  no  vertical  lift  to 
extend  the  trajectory.  A  profile  which  optimizes  the  crossrange  for  a  given  downrange 
is  needed. 

In  Reference  [10],  Spratlin  chose  to  fly  a  linear  bank  angle  with  velocity  profile 
for  entry  because  this  was  shown  to  nearly  maximize  crossrange  performance  in  [11] 
and  [12].  Wagner  [11]  used  several  different  optimization  techniques  to  show  that 
for  a  given  L/D,  the  bank  angle  profile  that  maximizes  crossrange  is  nearly  linear 
with  velocity.  This  can  be  seen  in  Figure  3.1  from  [11].  This  result  is  repeated  in 
Reference  [12]  with  a  varying  angle  of  attack  for  a  simulated  Shuttle  entry.  In  this 
study,  a  generalized  nonlinear  programming  technique  called  the  Hybrid  Nonlinear 
Programming  method  was  used  to  determine  the  optimal  bank  angle  profile  seen  in 
Figure  3.2.  Although  the  linear  bank  profile  is  not  optimal,  it  is  a  close  approximation 
of  the  optimal  profile  and  much  easier  to  implement.  An  additional  advantage  of  this 
profile  is  that  it  approaches  full  lift  up  at  the  end  of  the  trajectory  so  that  the 
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Figure  3.1.  Comparison  of  Optimal  Entry  Bank  Profiles 
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Bank  Angle,  de 


final  flight  path  angle  is  relatively  shallow,  and  this  is  one  of  the  desired  terminal 
conditions. 

Based  on  the  merits  of  the  linear  bank  profile,  it  was  chosen  to  be  used  to  demon¬ 
strate  vehicle  capabilities.  As  seen  in  Figure  3.3,  the  bank  angle  decreases  along  the 
profile  from  an  initial  value  at  entry  interface  to  zero  degrees  (full  lift- up)  at  the 
estimated  final  velocity. 

The  equation  which  defines  this  profile  is, 

<facmd  =  <t>des  ( - — )  (3.1) 

where, 

<i>cmd  is  the  commanded  bank  angle 
<j>des  is  the  initial  or  desired  bank  angle, 
v  is  the  current  inertial  velocity 
vei  is  the  entry  interface  velocity 
vj  is  a  previously  estimated  final  velocity 

Protection  is  provided  in  case  the  inertial  velocity  decreases  below  vj.  If  this  occurs, 
the  commanded  bank  angle  is  set  to  zero. 

The  initial  bank  angle  determines  the  amount  of  vertical  lift  which  will  be  used 
on  a  given  trajectory,  and  thus  defines  the  downrange  achievable  for  a  specific  linear 
bank  profile.  For  guidance  purposes,  the  profile  which  will  reach  a  target  downrange 
is  the  “desired”  profile,  so  the  initial  bank  angle  is  referred  to  as  the  desired  bank 
angle,  (fades- 
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3.3  General  Performance 

The  following  two  subsections  examine  the  impact  of  important  vehicle  and  environ¬ 
ment  factors  on  ranging  capability  and  final  dynamic  pressure.  The  vehicles  used  are 
generic  vehicles  defined  by  a  constant  L/D  and  Cb •  The  different  trajectories  flown 
are  obtained  by  varying  the  initial  bank  angle,  <j>dea. 

3.3.1  Ranging  Capability 

Crossrange  and  downrange  capability  results  are  presented  in  terms  of  a  maximum 
ranging  capability  footprint.  This  curve  demonstrates  the  furthest  point  in  a  given 
direction  that  a  vehicle  can  fly  using  the  near  optimal  linear  bank  profile.  The  max¬ 
imum  ranging  points  which  define  this  curve  were  obtained  by  flying  the  open-loop 
bank  profile  given  by  Equation  3.1  with  various  initial  bank  angles.  Figure  3.4  shows 
the  effect  of  L/D  on  ranging  capability.  Increased  L/D  has  a  large  impact  on  both 
crossrange  and  downrange  capabilities.  As  explained  in  Section  1.2,  crossrange  is  the 
more  limiting  factor  in  mission  flexibility.  For  a  vehicle  in  a  500  km  circular  orbit, 
there  is  approximately  1800  km  at  the  equator  between  adjacent  groundtracks  due 
to  planet  rotation.  This  requires  a  900  km  crossrange  capability  in  order  to  reach 
any  given  target  on  a  single  series  of  orbital  passes  and  ensure  a  deorbit  possibility 
within  a  reasonable  amount  of  time.  As  can  be  seen  in  Figure  3.4,  an  L/D  of  approx¬ 
imately  1.2  is  necessary  to  achieve  900  km  crossrange.  This  corresponds  to  an  L/D 
approaching  1.4  to  ensure  this  capability  in  the  presence  of  dispersions. 

The  effect  of  vehicle  ballistic  coefficient,  Cb,  on  ranging  capability  was  examined 
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Figure  3.4.  Maximum  Ranging  Capability:  Effect  of  L/D 
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next.  For  a  constant  L/D,  the  crossrange  capability  was  virtually  independent  of  bal¬ 
listic  coefficient,  and  the  downrange  performance  was  only  weakly  dependent  on  this 
parameter.  In  general,  vehicle  ranging  capability  is  not  effected  by  Cb  as  supported 
by  Figure  3.5. 

The  variation  in  ranging  capability  due  to  the  coefficient  of  lift  is  the  same  as 
the  variation  due  to  L/D  because  changing  Ci  simply  effects  L/D.  The  variation  in 
ranging  capability  due  to  the  coefficient  of  drag  is  almost  the  same  as  the  variation 
due  to  L/D.  This  is  true,  despite  the  fact  that  changing  Co  also  changes  Cb,  because 
ranging  capability  is  essentially  independent  of  Cb- 

The  only  initial  condition  for  which  ranging  performance  is  found  to  be  sensitive 
is  the  entry  flight  path  angle,  7*,-.  Small  changes  in  initial  flight  path  angle  greatly 
vary  downrange  performance,  and  also  have  a  significant,  although  smaller,  impact 
on  crossrange  capability.  As  seen  in  Figure  3.6,  the  downrange  performance  can  be 
greatly  increased  by  entering  at  a  shallower  flight  path  angle.  However,  this  increase  in 
downrange  performance,  which  only  slightly  enhances  guidance  capability,  comes  with 
an  increased  risk  of  the  vehicle  skipping  back  out  of  the  atmosphere.  The  maximum 
downrange  (full  lift-up)  case  with  the  shallow  7e,  came  very  close  to  skipping  out, 
and  any  shallower  flight  path  angle  would  probably  skip.  The  crossrange  is  increased, 
but  the  risk  associated  with  skipping  and  widely  missing  the  target  must  be  balanced 
with  this. 

The  steeper  7 ei  reduced  both  the  downrange  and  crossrange  maneuver  capacity, 
but  the  change  in  ranging  capability  is  not  as  large  as  it  is  for  the  shallower  7 ej. 
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Figure  3.5.  Maximum  Ranging  Capability:  Effect  of  Ballistic  Coefficient 
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Figure  3.6.  Maximum  Ranging  Capability:  Effect  of  Entry  Flight  Path  Angle 
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However,  the  difference  is  still  s'gnificant  as  the  reduction  is  approximately  equivalent 
to  a  reduction  in  L/D  from  1.0  to  0.9.  The  large  variation  in  ranging  capability  for 
such  small  differences  in  7e,-  indicates  that  it  will  probably  be  important  to  hit  the 
target  value  because  there  is  a  narrow  entry  corridor.  The  target  entry  flight  path 
angle  should  be  chosen  as  a  compromise  between  maximizing  ranging  capability  and 
minimizing  the  risk  of  skipping. 

The  two  major  environmental  concerns  during  entry  are  atmospheric  density  and 
wind  dispersions.  Winds  do  not  significantly  affect  the  ranging  performance  as  they 
might  be  expected  to.  For  example,  a  headwind  increases  the  lift  of  the  vehicle  so 
that  it  flies  for  a  greater  length  of  time,  and  this  partially  counters  the  effect  reduced 
ground  speed  has  on  range.  Crossrange  is  affected  more  than  downrange  because  the 
vehicle  crossrange  velocity  is  slower  relative  to  the  wind  than  the  downrange  velocity. 
The  decrease  in  crossrange  caused  by  an  opposing  wind  is  noticeable,  but  not  a  major 
ranging  factor. 

Atmospheres  which  have  variable  densities  that  continually  shift  about  the  nomi¬ 
nal  density  have  a  relatively  small  effect  on  ranging  capability  because  the  integrated 
effect  is  about  the  same  as  if  the  density  had  remained  nominal.  (While  vehicle  per¬ 
formance  is  not  greatly  affected,  this  type  of  atmosphere  can  profoundly  influence 
guidance  performance,  as  will  be  examined  in  Chapter  5.)  It  is  large  density  biases 
which  remain  throughout  the  trajectory  that  can  have  a  significant  effect  of  ranging 
performance.  This  effect  can  be  seen  in  Figure  3.7.  An  atmosphere  which  is  consis¬ 
tently  thinner  than  expected  increases  vehicle  downrange  capability  because  there  is 
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Figure  3.7.  Maximum  Ranging  Capability:  Effect  of  Constant  Density  Biases 
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less  resistance,  but  at  the  same  time,  this  atmosphere  decreases  crossrange  capability 
because  there  is  less  atmosphere  to  turn  the  vehicle.  A  thick  atmosphere  decreases 
downrange  performance  and  increases  crossrange  capability  for  the  opposite  reasons. 
It  is  interesting  to  note  that  a  headwind  acts  much  like  a  thick  atmosphere  and  a 
tailwind  acts  like  a  thin  atmosphere  in  regards  to  ranging  performance,  the  difference 
being  that  the  change  due  to  winds  is  much  smaller  than  the  change  due  to  density 
biases  for  the  assumed  dispersions. 

3.3.2  Final  Dynamic  Pressure 

The  influence  of  vehicle  design  and  environmental  factors  on  final  dynamic  pressure, 
qf ,  is  examined  in  this  section.  But  first,  the  nature  of  the  dynamic  pressure  variation 
during  entry  is  presented. 

The  dynamic  pressure  history  throughout  entry  was  found  to  be  analogous  to  a 
damped  sinusoid  as  can  be  seen  in  an  example  in  Figure  3.8.  This  is  not  unexpected  as 
the  dynamics  of  a  constant  L/D  entry  trajectory  have  been  likened  to  the  response  of  a 
simple  second-order  spring-mass-damper  system  in  [13].  During  the  trade-off  between 
potential  and  kinetic  energy  which  occurs  throughout  entry,  drag  is  a  non-conservative 
force  which  provides  damping  to  the  system.  Hence,  the  dynamic  pressure  appears 
to  be  driving  to  a  fmai  "steady-state”  value  with  long-period  (phugoid)  oscillations 
about  that  value.  The  final  dynamic  pressure  varies  widely  for  a  given  vehicle  design  or 
dispersion  because  the  magnitude  of  qj  is  dependent  not  only  upon  the  “steady-state” 
value,  but  upon  the  point  in  the  oscillation  about  that  value  that  the  entry  maneuver 
terminates.  The  peaks  and  valleys  in  the  typical  trajectory  shown  in  Figure  3.8 
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DYNAMIC  PRESSURE  (N/SQ  Ml 


range  between  1150  and  800  N/sq  m  while  the  “steady- state”  dynamic  pressure  is 
approximately  1000  N/sq  m. 

It  is  clear  that  the  final  dynamic  pressure  for  a  single  trajectory  is  not  sufficient 
to  define  the  qj  qualities  associated  with  a  given  vehicle  design  or  dispersion.  For 
this  reason,  the  range  of  final  dynamic  pressures  for  a  wide  selection  of  trajectories 
is  necessary  to  characterize  dynamic  pressure  performance.  The  various  trajectories 
used  to  define  the  maximum  ranging  capabilities  in  the  previous  section  conveniently 
provide  a  database  of  final  dynamic  pressures  which  is  used  in  this  analysis.  The 
range  of  final  dynamic  pressures  for  a  given  parameter  represents  the  magnitude  of 
the  oscillations  about  the  “steady-state”  value,  while  the  median  value  in  the  range 
is  a  rough  indicator  of  the  “steady-state”  final  dynamic  pressure. 

The  first  vehicle  characteristic  analyzed  for  dynamic  pressure  performance  was 
L/D.  Increasing  L/D  was  found  to  decrease  final  dynamic  pressure  as  indicated  in 
Figure  3.9.  A  higher  L/D  allows  the  vehicle  to  fly  further  through  the  atmosphere,  and 
this  means  more  atmosphere  to  decelerate  the  vehicle.  This  implies  that  the  vehicle 
design  L/D  should  be  as  large  as  feasible  if  final  dynamic  pressure  is  a  concern. 

Decreasing  ballistic  coefficient  was  found  to  also  decrease  final  dynamic  pressure 
as  seen  in  Figure  3.10.  The  mean  value,  as  well  as  the  size  of  the  range  in  values, 
decreased  with  Cb •  In  fact,  for  any  given  trajectory  defined  by  its  bank  profile, 
the  final  dynamic  pressure  was  found  to  be  almost  linearly  dependent  upon  Cb-  This 
result  is  demonstrated  for  a  near-maximum  crossrange  trajectory  which  used  an  initial 
bank  angle  of  65°  in  Figure  3.11  .  This  implies  that  the  entry  vehicle  design  should 
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Figure  3.9.  Final  Dynamic  Pressure  Range:  Effect  of  L/D 


62 


Figure  3.10.  Final  Dynamic  Pressure  Range:  Effect  of  Ballistic  Coefficient 
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FINAL  DYNAMIC  PRESSURE 
(N/*q  m) 


have  as  low  as  ballistic  coefficient  as  possible  when  final  dynamic  pressure  is  a  factor. 

Next,  the  environmental  effects  on  final  dynamic  pressure  are  analyzed.  Winds 
were  found  to  be  a  minor  factor  influencing  final  dynamic  pressure.  Regardless  of  the 
winds,  the  vehicle  generally  reached  the  target  with  approximately  the  same  relative 
velocity,  and  hence  the  same  qj.  In  fact,  the  primary  environment  characteristic 
affecting  final  dynamic  pressure  is  density.  An  atmosphere  with  a  constant  thin  bias 
does  not  permit  generation  of  enough  drag  to  allow  the  vehicle  to  sufficiently  slow 
down,  causing  higher  final  dynamic  pressures.  A  thick  atmosphere  has  the  opposite 
effect.  It  allows  the  vehicle  to  decelerate  quickly  and  helps  reduce  final  dynamic 
pressure.  The  dramatic  effect  variations  in  atmospheric  density  have  on  final  dynamic 
pressure  is  evident  in  Figure  3.12.  This  is  a  concern  because  there  is  a  relatively  large 
amount  of  uncertainty  regarding  the  density  that  will  be  encountered. 

In  the  last  portion  of  this  study,  the  effect  that  dispersions  in  the  entry  flight 
path  angle  had  on  final  dynamic  pressure  was  examined.  The  results  are  presented  in 
Figure  3.13.  The  steeper  flight  path  angle  caused  a  wider  range  in  the  final  dynamic 
pressure,  but  the  mean  value  remained  the  same  regardless  of  the  variation.  Some 
caution  must  be  used  in  interpreting  this  because  these  results  are  specific  for  a 
single  7e,.  Other  entry  flight  path  angles  might  be  more  sensitive,  in  terms  of  qj,  to 
dispersions. 
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Figure  3.13.  Final  Dynamic  Pressure  Range:  Effect  of  Entry  Flight  Path  Angle 
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3.4  Vehicle  M  Performance 


The  design  of  entry  vehicles  is  not  driven  by  guidance  considerations  alone.  Once  a 
vehicle  has  been  designed,  a  guidance  must  make  do  with  the  vehicle  capabilities  and 
limitations.  So,  given  a  particular  vehicle,  it  is  important  to  determine  how  the  per¬ 
formance  can  be  varied  and  what  the  maximum  possible  performance  is.  The  vehicle 
that  was  used  to  demonstrate  this  analysis  was  a  biconic  lifting  body,  designated 
Vehicle  M,  baselined  at  NASA  Johnson  Space  Center  for  the  Mars  Rover  Sample 
Return  program  [1].  Vehicle  M  has  a  maximum  L/D  near  1.1,  and  at  the  trim  angle 
of  attack  of  26.04°,  Vehicle  M  has  a  1.0  L/D.  The  lift  coefficient,  drag  coefficient, 
and  L/D  curves  which  define  this  vehicle  can  be  seen  in  Figures  3.14,  3.15,  and  3.16. 

In  addition,  the  mass  of  Vehicle  M  is  4082  kg,  and  the  surface  reference  area  is 
14.19  sq  m.  This  reference  area  combines  with  the  mass  and  the  coefficient  of  drag 
at  the  trim  angle  of  attack  to  give  a  ballistic  coefficient  of  397  kg/sq  m. 

For  Vehicle  M,  the  trim  angle  of  attack  is  very  sensitive  to  the  location  of  the 
center  of  gravity.  Near  30°,  a  1%  shift  in  center  of  gravity  location  along  the  length 
of  the  vehicle  will  cause  a  2.5°  shift  in  trim  angle  of  attack  [14].  Therefore,  it  is  a 
relatively  simple  matter  to  change  a*.  Since  changing  at  varies  the  effective  ballistic 
coefficient  as  well  as  changing  the  lift-to-drag  ratio,  the  vehicle  performance  can  be 
modified  by  small  changes  in  the  placement  of  the  center  of  gravity. 

The  previous  section  showed  that  ranging  capability  is  basically  defined  by  the 
L/D  at  which  the  vehicle  flies — the  higher  the  L/D,  the  greater  the  ranging  capability. 
To  maximize  the  ranging  capability  of  Vehicle  M,  the  vehicle  should  be  flown  at  the 
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Figure  3.14.  Vehicle  M  Lift  Coefficient  Curve 
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Figure  3.15.  Vehicle  M  Drag  Coefficient  Curve 
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Figure  3.16.  Vehicle  M  L/D  Curve 
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trim  angle  of  attack  which  corresponds  to  maximum  L/D.  It  is  also  noted  that  if 
one  is  not  flying  at  the  maximum  L/D,  it  is  possible  to  fly  the  same  L/D,  and  get 
approximately  the  same  ranging  performance,  at  two  different  angles  of  attack.  This 
is  illustrated  for  the  trim  angles  of  attack  of  12°  and  26.04°  in  Figure  3.17.  The 
difference  in  ranging  performance  is  just  the  small  change  caused  by  variation  of 
ballistic  coefficient  that  was  noted  in  Section  3.3.1. 

However,  ranging  capability  is  not  the  only  performance  requirement.  Final  dy¬ 
namic  pressure  performance  must  also  be  considered  when  determining  the  trim  angle 
of  attack  to  be  flown.  The  change  in  final  dynamic  pressure  due  to  various  trim  angles 
of  attack  for  Vehicle  M  is  shown  in  Figure  3.18.  It  can  be  seen  that  the  higher  angles 
of  attack  result  in  lower  final  dynamic  pressures.  This  is  because  the  higher  angles 
of  attack  increase  the  coefficient  of  drag  which  has  the  same  effect  on  final  dynamic 
pressure  as  reducing  the  ballistic  coefficient.  The  reduction  in  final  dynamic  pressure 
(and  final  velocity)  obtained  by  flying  at  the  higher  angles  of  attack  is  the  reason 
that  entry  vehicles  generally  fly  on  the  back  side  of  the  L/D  curve.  The  back  side  of 
the  L/D  curve  is  where  increasing  angle  of  attack  results  in  lower  lift-to-drag  ratios 
as  seen  in  Figure  3.16.  The  Space  Shuttle  flies  on  the  back  side  of  its  L/D  curve  until 
late  in  the  trajectory  when  it  transitions  from  a  decelerating  entry  vehicle  to  a  flying 
glider. 

Besides  final  dynamic  pressure  considerations,  there  are  several  other  reasons  why 
it  is  desirable  to  fly  at  higher  angles  of  attack.  A  higher  angle  of  attack  means  a 
further  aft  center  of  gravity  (eg)  location.  Reference  [14]  states  that  a  further  aft 
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CROSSRANGE 


Figure  3.17.  Vehicle  M:  Ranging  Performance  Comparison  of  Two  Angles  of  Attack 
—  Front  and  Back  Side  of  L/D  Curve 


73 


74 


eg  is  advantageous  because  1)  during  launch,  the  eg  is  lower  and  closer  to  the  base 
which  provides  better  booster  performance  due  to  reduced  booster  loads  and  2)  it 
is  much  easier  to  package  components  in  the  large  end  of  the  vehicle  cone.  While 
the  first  advantage  is  specific  to  one  candidate  booster,  the  Centaur  upper  stage,  the 
second  is  true  of  most  any  biconic  vehicle.  Additionally,  a  higher  angle  of  attack  is 
statically  more  stable  which  means  that  there  is  greater  restoring  force  if  the  vehicle 
is  perturbed  from  the  nominal  angle  of  attack.  A  final  reason  why  a  higher  angle 
of  attack  is  preferable  is  that  it  decreases  the  ballistic  lift  coefficient,  m/CxS,  which 
means  that  the  vehicle  is  more  maneuverable  and  can  effect  faster  turn  rates. 

While  all  of  these  factors  favor  a  high  angle  of  attack,  it  must  be  realized  that 
on  the  back  side  of  the  L/D  curve,  L/D  and  ranging  capability  decrease  with  greater 
angles  of  attack.  Because  of  this  trade-off,  Vehicle  M  should  be  flown  only  at  a 
high  enough  angle  of  attack  to  sufficiently  reduce  q/.  Further  increases  in  at  are 
not  justified  in  view  of  decreasing  ranging  performance.  This  decrease  is  quantified 
for  Vehicle  M  in  Figure  3.19.  For  the  given  constraints,  an  angle  of  attack  of  33°  is 
necessary  to  reduce  the  final  dynamic  pressure  to  within  the  general  limits.  As  can  be 
seen,  this  causes  a  significant  reduction  in  ranging  performance,  illustrating  the  trade¬ 
off  for  fixed  vehicle  aerodynamic  design  between  maximizing  ranging  capabilities  and 
decreasing  final  dynamic  pressure. 
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Figure  3.19.  Vehicle  M:  Effect  of  Angle  of  Attack  on  Maximum  Ranging  Capability 


3.5  Divert  Capability 

During  entry,  there  will  undoubtedly  be  inaccuracy  in  the  vehicle  navigation  position 
estimate  and  in  the  knowledge  of  the  target  location  due  to  significant  errors  in 
knowledge  of  the  orientation  of  Mars  in  inertial  space  and  in  the  ability  to  locate 
remote  images  on  the  globe.  This  will  cause  the  guidance  to  aim  for  a  point  which 
is  not  actually  the  true  target.  In  order  to  reduce  these  errors,  a  beacon  may  be 
placed  at  the  target.  Navigation  updates  provided  by  this  beacon  will  help  determine 
the  location  of  the  vehicle  relative  to  the  landing  site,  and  in  doing  so  will  “move” 
the  target  relative  to  the  vehicle.  The  vehicle  must  still  be  able  to  reach  the  target 
even  though  it  was  not  initially  aiming  for  it.  This  matter  is  complicated  because  this 
navigation  aid  only  becomes  visible  late  in  the  entry  trajectory.  Therefore,  the  ability 
of  a  vehicle  to  divert  to  the  new  target  location  late  in  the  trajectory  is  important. 

In  order  to  investigate  divert  capability,  an  entry  was  flown  to  a  typical  target 
(DR  3500  km,  CR  200  km).  For  this  trajectory,  the  landing  site  became  visible  at 
a  range  of  468  kilometers.  It  was  assumed  for  this  study  that  the  navigation  beacon 
which  provides  improved  knowledge  of  the  landing  site  location  also  becomes  visible 
at  this  point.  From  this  possible  divert  point  in  the  nominal  trajectory,  the  open-loop 
ranging  capabilities  of  Vehicle  M  (at  =  33°,  L/D  =  0.85)  were  analyzed.  In  addition, 
the  divert  capability  was  tested  at  t  vo  later  points  in  the  trajectory — 278  and  103  km 
from  the  target.  The  inertial  velocity  at  each  of  these  three  points  was  approximately 
1400,  1000,  and  500  m/s  greater  than  the  nominal  terminal  velocity.  The  velocity 
prior  to  nominal  termination  is  given  because  it  is  a  more  universal  indicator  of  the 
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location  of  the  vehicle  in  the  trajectory.  The  results  are  shown  in  Figure  3.20. 

The  footprint  shrinks  rapidly  as  the  vehicle  nears  the  target.  Immediately  after 
signal  acquisition  (1400  m/s  prior  to  ending),  under  ideal  conditions  the  vehicle  could 
reach  any  new  location  within  a  45  km  radius  of  the  original  target.  Approximately 
190  km  and  400  m/s  later,  the  footprint  has  shrunk  to  20  km.  Another  175  km  and 
500  m/s  later,  at  a  range  of  103  km,  the  circle  radius  is  only  7  km. 

The  larger  divert  bank  angles  cause  the  entry  vehicle  to  dump  a  significant  amount 
of  vertical  lift  in  order  to  reach  a  closer  target.  Since  a  loss  in  vertical  lift  initially 
causes  a  drop  in  altitude  and  an  increase  in  velocity,  the  final  dynamic  pressure  can 
be  significantly  increased  by  attempts  to  divert  to  short  targets  late  in  the  trajectory. 
The  increase  in  final  dynamic  pressure  is  not  as  great  when  the  divert  is  earlier  in  the 
trajectory  because  the  increased  drag  due  to  the  higher  velocity  and  density  has  time 
to  act  on  the  vehicle  and  reduce  velocity.  For  the  trajectory  tested,  the  vehicle  could 
divert  anywhere  within  its  ranging  capability  at  signal  acquisition  without  exceeding 
the  final  dynamic  pressure  limit.  For  the  later  diverts,  some  shorter  targets  are  not 
reachable  because  of  the  final  dynamic  pressure  limit.  If  divert  capability  is  defined 
as  the  largest  possible  circle  centered  at  the  target  which  encompasses  divert  landing 
sites,  the  dynamic  pressure  limit  does  not  really  effect  the  divert  capability.  This  is 
true  because  the  size  of  this  circle  is  bounded  in  the  positive  downrange  direction  by 
a  limited  ability  to  fly  further. 

Without  any  navigation  updates  from  an  orbiting  vehicle  or  other  ground  bea¬ 
cons  during  entry,  the  3 <r  navigation  position  error  relative  to  the  landing  site  at 
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Figure  3.20.  Divert  Capability 


79 


parachute  deployment  is  expected  to  be  less  than  5  km  [15].  Provided  that  the  need 
to  divert  is  recognized  soon  after  acquisition  of  signal,  the  divert  capability  for  this 
low  performance  vehicle  is  well  within  expected  navigation  errors. 
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Chapter  4 

Guidance  Design 


The  guidance  algorithm  developed  here  uses  bank  angle  control  in  order  to  meet 
downrange,  crossrange,  and  final  dynamic  pressure  constraints.  This  chapter  exam¬ 
ines  previous  work  relevant  to  entry  guidance  in  Section  4.1  and  then  presents  the 
guidance  and  its  motivation  in  Section  4.2.  The  derivation  of  each  portion  of  the 
guidance  algorithm  follows  in  the  remainder  of  the  chapter  with  the  downrange  con¬ 
trol  explained  in  Section  4.3,  the  crossrange  control  detailed  in  Section  4.4,  and  the 
dynamic  pressure  control  examined  in  Section  4.5. 

4.1  Previous  Work 

The  majority  of  precision  landing  guidance  work  has  been  concentrated  on  Earth 
entry,  and  relatively  little  guidance  work  has  been  done  regarding  precision  entry  at 
Mars.  Although  there  are  significant  differences,  such  as  a  thinner,  more  uncertain 
Mars  atmosphere  and  greater  navigation  errors,  the  two  problems  are  similar  enough 
for  much  of  the  Earth  work  to  be  applicable  to  Mars.  Earth  entry  guidance  designs 
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which  are  of  interest  include  those  for  Apollo  [16],  the  Space  Shuttle  Orbiter  [17],  and 
the  Entry  Research  Vehicle  (ERV)  [10].  Also  of  interest  is  an  aerobraking  guidance 
algorithm  developed  for  the  Aeroassist  Flight  Experiment  (AFE)  [18]. 

The  Apollo  entry  guidance  was  designed  to  obtain  landing  accuracy  of  15  nautical 
miles  using  a  very  low  L/D  (0.3)  vehicle  [16].  This  was  achieved  using  one  control 
variable,  bank  angle,  to  maintain  a  reference  drag  acceleration  profile  trajectory. 
This  reference  drag  profile  was  analytically  determined  prior  to  flight  to  meet  the 
downrange,  heating,  and  g-load  requirements.  The  crossrange  requirements  were  met 
by  varying  the  sign  of  the  bank  angle  based  on  the  current  crossrange  error.  The 
low  L/D  of  the  vehicle,  while  providing  relatively  little  control  over  the  trajectory, 
also  reduced  the  length  of  the  trajectory.  The  shorter  trajectory  reduced  the  loss  in 
accuracy  due  to  poor  control  authority  by  subjecting  the  vehicle  to  fewer  atmospheric 
dispersions. 

The  Shuttle  entry  guidance  is  designed  to  guide  a  mid  L/D  vehicle  (L/D  «  1.2) 
from  a  low-Earth  orbit  to  a  runway  landing.  Although  the  Shuttle  guidance  has  a 
much  more  aerodynamic  vehicle  to  work  with  than  Apollo,  the  accuracy  requirements 
are  much  greater.  The  Shuttle  guidance  uses  three  control  variables,  bank  angle,  angle 
of  attack,  and  a  speed  brake,  to  fly  reference  drag  and  altitude  rate  profiles.  These 
reference  profiles  are  analytically  determined  to  attain  the  desired  downrange  and  to 
avoid  flight  conditions  which  approach  the  heating  and  g-load  limits.  Assumptions 
concerning  the  vehicle  performance  and  environment  are  used  in  determining  these 
profiles  prior  to  flight.  The  crossrange  is  controlled  by  bank  reversals  similar  to 
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Apollo,  but  the  reversals  are  based  on  the  difference  between  the  current  heading  and 
the  heading  to  the  target  as  opposed  to  the  crossrange  error. 

Although  not  an  entry  guidance,  Higgins’  candidate  aerobraking  guidance  for  the 
Aeroassist  Flight  Experiment  in  [18]  utilizes  many  concepts  applicable  to  entry.  The 
AFE  vehicle  enters  the  Earth’s  atmosphere  from  a  high  energy  orbit  comparable  to  a 
geosynchronous  return  orbit  and  skims  through  the  upper  portions  of  the  atmosphere 
to  reduce  its  energy.  It  then  exits  the  atmosphere  into  a  low  Earth  orbit.  The 
aerobraking  guidance  uses  only  one  control  variable,  bank  angle,  to  target  an  apogee 
and  plane  for  this  orbit.  A  target  orbital  apogee  is  similar  to  a  target  downrange, 
while  plane  change  capability  is  analogous  to  crossrange  capability.  This  is  because 
the  desired  apogee  is  dependent  upon  the  magnitude  of  the  in-plane  lift  just  like 
downrange,  and  the  change  in  orbital  plane  is  effected  by  the  lateral  lift  just  like 
crossrange.  This  guidance  uses  the  bank  angle  to  control  the  magnitude  of  the  in¬ 
plane  lift  to  target  to  a  desired  apogee.  Plane  errors  are  corrected  using  any  residual 
out-of-plane  lift,  and  the  direction  of  the  out-of-plane  lift  is  controlled  by  the  sign  of 
the  bank  angle.  In  this  manner,  the  one  control  variable,  bank  angle,  is  used  to  affect 
two  conditions. 

In  Higgins’  guidance,  the  bank  angle  necessary  to  reach  a  target  apogee  is  deter¬ 
mined  by  a  numeric  predictor-corrector  algorithm.  A  predictor-corrector  algorithm 
computes  control  corrections  based  on  predictions  of  the  terminal  conditions  obtained 
by  flying  an  assumed  control  profile.  The  use  of  a  predictor-corrector  implies  a  ref¬ 
erence  control  profile  (in  the  case  of  AFE,  a  constant  bank  angle  profile),  as  opposed 
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to  a  reference  trajectory.  Reference  trajectories,  such  as  were  used  in  the  Apollo  and 
Shuttle  guidances,  are  excellent  at  controlling  the  vehicle  energy  state  and  reaching 
a  target  when  there  are  few  unknowns.  However,  this  form  of  guidance  can  have 
difficulties  if  the  vehicle  deviates  too  far  from  the  reference  trajectory.  This  could 
happen  if  there  are  large  variations  from  the  expected  conditions  during  entry,  as  is 
likely  at  Mars. 

As  implemented  in  the  proposed  AFE  guidance,  the  predictor-corrector,  using  a 
reference  bank  profile,  essentially  recomputes  a  new  reference  trajectory  each  guid¬ 
ance  cycle,  and  thus  is  able  to  better  adapt  to  wide  variations  in  atmospheric  and 
vehicle  conditions.  A  drawback  of  the  reference  control  profile  is  that  it  guides  the 
vehicle  based  on  the  predicted  final  state  and  does  not  concern  itself  with  the  interme¬ 
diate  states.  The  state  of  the  vehicle  at  intermediate  steps  can  be  just  as  important 
as  the  terminal  conditions,  as  there  are  continuous  constraints,  such  as  g-load  and 
heating  rate,  which  the  guidance  can  violate  while  still  achieving  the  desired  terminal 
conditions. 

The  prediction  algorithm  used  in  the  guidance  could  have  been  either  analytic  or 
numeric.  An  analytic  predictor  uses  closed-form  solutions  of  the  equations  of  motion 
to  predict  the  final  state.  However,  in  order  to  obtain  the  closed-form  solutions,  some 
simplifying  assumptions  must  be  made.  This  reduces  the  accuracy  of  the  analytic  pre¬ 
dictions.  The  numeric  predictor,  which  was  chosen  by  Higgins  for  AFE,  numerically 
integrates  the  equations  of  motion  forward  in  time  using  assumed  atmosphere  and  ve¬ 
hicle  models  and  obtains  highly  accurate  predictions.  The  cost  of  this  accuracy  is  an 
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increased  computational  load.  A  numeric  predictor- corrector  can  be  very  inefficient 
and  computer-resource  intensive  if  improperly  designed.  The  numeric  algorithm  is 
now  viable  because  computational  abilities  have  greatly  improved  in  the  last  decade, 
and  more  importantly,  the  basic  algorithm  was  refined  by  Spratlin. 

Spratlin  recognized  the  capabilities  of  the  numeric  predictor-corrector  type  algo¬ 
rithm  and  applied  this  concept  to  entry  with  the  Entry  Research  Vehicle  in  Ref¬ 
erence  [10].  The  ERV  was  proposed  to  demonstrate  synergetic  plane  change  and 
precision  entry  capabilities  with  a  high  L/D  (%  2.0)  vehicle  at  Earth.  The  ERV 
guidance  is  based  upon  linear  bank  angle  with  velocity  and  constant  angle  of  attack 
control  profiles.  A  numeric  predictor-corrector  determines  the  change  in  bank  angle 
and  angle  of  attack  necessary  to  reach  the  desired  crossrange  and  downrange.  The 
commands  are  computed  by  solving  a  pair  of  simultaneous  equations  involving  bank 
angle,  angle  of  attack,  and  the  sensitivity  of  these  controls  to  the  currently  predicted 
downrange  and  crossrange  errors.  The  heating  rate  is  controlled  explicitly  during  the 
short  periods  when  heating  is  a  problem.  No  attempt  is  made  to  control  the  energy 
level  of  the  entry  vehicle  at  the  target. 

A  reference  profile,  if  not  closely  followed,  quickly  diverges  from  the  desired  tra¬ 
jectory,  so  errors  must  be  continuously  fed  back  to  recompute  the  commanded  bank 
angle  in  an  analytic  reference  trajectory  guidance.  Spratlin  reasoned  that  the  nu¬ 
meric  predictor-corrector  in  the  ERV  guidance  need  not  be  called  as  frequently  be¬ 
cause  the  trajectory  solution  (the  reference  control  profile)  does  not  change  as  rapidly. 
The  commanded  bank  angle  is  still  determined  every  guidance  cycle  from  the  cur- 
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rent  reference  bank  profile,  but  a  new  reference  control  profile  is  recomputed  by  the 
predictor-corrector  at  a  much  slower  rate.  This  decreases  the  computational  load. 
An  additional  improvement  in  the  algorithm  is  the  steps  taken  towards  optimizing 
the  guidance  rate.  The  rate  could  be  varied  based  on  aerodynamic  accelerations  or 
time-to-go  to  produce  a  nearly  constant  load  on  the  central  processing  unit.  This 
allows  the  guidance  scheme  to  be  run  at  a  slow  rate  early  in  the  trajectory  when 
the  trajectory  to  be  predicted  is  long  and  there  is  plenty  of  control  margin.  Later, 
when  the  time-to-go  is  short  and  control  authority  is  small,  thereby  making  frequent 
guidance  updates  necessary,  the  guidance  should  be  run  at  a  faster  rate. 

Previous  Mars  entry  guidance  studies  have  largely  been  devoted  to  low  L/D  entries 
such  as  those  used  by  the  Viking  landers  [19].  However,  in  Reference  [5],  Tigges  has 
recently  developed  a  precision  entry  scheme  for  a  mid  L/D  vehicle  at  Mars.  The 
guidance  analytically  predicts  the  equilibrium  glide  profile  which  gives  the  desired 
downrange.  The  crossrange  is  controlled  by  bank  reversals  based  on  the  current 
crossrange  error.  The  guidance  consists  of  two  phases.  The  first  phase  uses  bank  angle 
control  only  to  meet  the  ranging  constraints.  Late  in  the  trajectory,  if  the  dynamic 
pressure  limit  is  predicted  to  be  exceeded,  then  the  second  phase  is  implemented.  In 
the  second  phase,  angle  of  attack  modulation  is  used  in  addition  to  bank  angle  control 
to  meet  the  ranging  and  dynamic  pressure  constraints. 


86 


4.2  Guidance  Scheme 


Precision  landing  requires  robust  control  over  the  vehicle  trajectory.  As  demonstrated 
in  [5],  the  necessary  control  can  be  achieved  at  Mars  by  varying  both  the  vehicle  bank 
angle  and  angle  of  attack.  However,  because  the  ability  to  modify  angle  of  attack 
adds  complexity  to  an  entry  vehicle  (as  explained  in  Section  1.3),  it  is  desirable  to 
achieve  precision  entry  using  a  vehicle  with  bank  angle  control  only.  Thus,  this 
study  has  chosen  to  simplify  the  vehicle  by  reducing  the  control  variables.  The  one 
remaining  control  variable  has  two  components — the  magnitude  and  sign  of  the  bank 
angle — which  must  be  used  to  reach  the  target  and  satisfy  the  constraints. 

The  constraints  which  must  be  met  have  been  listed  in  Table  1.1.  However,  not 
all  of  the  constraints  must  be  actively  controlled.  For  the  entry  initial  conditions 
chosen  (i.e.  deorbit  from  a  low  circular  orbit),  g-ioading  and  heating  rate  are  not  a 
major  concern  because  the  g-loads  experienced  during  these  entries  are  less  than  those 
expected  upon  launch  at  Earth,  and  the  heating  encountered  is  relatively  benign  and 
well  within  the  capabilities  of  state  of  the  art  thermal  protection  systems.  Active 
control  of  the  final  flight  path  angle  was  initially  considered,  but  was  later  found 
unnecessary  because  shallow  final  flight  path  angles  naturally  resulted  from  the  chosen 
guidance  design. 

The  major  constraints  which  the  guidance  must  meet  are  final  range  from  the 
target  and  final  dynamic  pressure.  The  final  range  error  can  be  expressed  as  cross¬ 
range  and  downrange  components.  The  two  ranging  components  are  assumed  to  be 
independently  controllable  in  the  guidance  design,  with  the  root-sum-squared  of  the 
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two  errors  at  guidance  termination  indicating  ranging  performance.  This  results  in  a 
total  of  three  constraints  which  the  guidance  must  actively  meet. 

The  guidance  must  not  violate  any  of  the  listed  constraints  despite  dispersions  in 
the  atmosphere  and  vehicle  performance.  The  major  dispersions  which  the  guidance 
must  be  prepared  to  handle  are  off-nominal  atmospheric  densities,  winds,  and  erro¬ 
neous  predictions  in  vehicle  aerodynamics.  The  guidance  must  also  maximize  cross¬ 
range  capability  and  conserve  total  ranging  capability.  As  explained  in  Section  1.2, 
the  crossrange  capability  is  important  to  increase  mission  flexibility.  The  ability  to 
conserve  total  ranging  capability  is  also  important  because  it  is  an  indication  of  the 
ability  of  the  guidance  to  handle  dispersions  later  in  the  trajectory. 

Prior  to  determining  the  general  guidance  scheme,  it  is  necessary  to  know  how 
the  one  control,  bank  angle,  influences  the  constraints.  The  bank  angle  determines 
the  orientation  of  the  lift  vector  in  the  plane  perpendicular  to  the  relative  velocity, 
and  varying  this  angle  apportions  the  lift  between  the  in-plane  (or  vertical)  and  out- 
of-plane  (or  lateral)  directions.  The  vertical  lift  affects  both  the  downrange  and 
final  dynamic  pressure.  While  the  downrange  distance  is  directly  influenced  by  the 
vertical  lift  throughout  the  trajectory,  the  final  dynamic  pressure  is  not  so  obviously  a 
result  of  the  vertical  lift.  As  discussed  in  Section  3.3.2,  the  dynamic  pressure  history 
during  entry  approximates  a  second-order  system  with  a  phugoid  mode,  and  the  final 
dynamic  pressure  is  determined  by  the  interaction  of  a  number  of  factors,  including 
bank  profile  (i.e.  the  vertical  lift).  The  lateral  lift  affects  the  crossrange  distance  to  be 
flown.  Greater  lateral  lift,  given  the  same  vertical  lift,  results  in  greater  crossrange. 
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However,  vertical  lift  is  also  an  important  factor  in  crossrange  because  it  determines 
how  long  the  lateral  lift  has  to  act. 

Since  there  are  three  constraints  and  only  one  control  variable,  trade-offs  must 
be  made.  Design  of  the  guidance  begins  with  the  downrange  control.  After  control 
of  the  downrange  constraint  is  established,  control  of  the  other  two  constraints  will 

follow. 

Downrange  control  is  accomplished  using  a  numeric  predictor-corrector  algorithm 
which  controls  the  magnitude  of  the  current  bank  angle  based  on  the  error  in  the 
final  downrange  obtained  by  flying  an  assumed  bank  history.  The  reference  control 
profile  and  predictor-corrector  algorithm  were  chosen  over  a  reference  trajectory  guid¬ 
ance  scheme  based  on  the  uncertainty  in  the  Mars  environment  and  the  need  for  the 
adaptability  that  the  predictor-corrector  provides.  The  prediction  was  chosen  to  be 
performed  numerically  because  of  the  increased  accuracy  over  analytic  predictions. 
Although  not  done  here,  the  predictor-corrector  execution  rate  and  the  prediction 
step  size  can  both  be  optimized  to  reduce  the  computational  load.  With  these  up¬ 
grades,  the  computational  requirements  for  the  numeric  predictor-corrector  guidance 
are  not  excessive,  and  the  scheme  is  a  viable  option  for  a  Mars  mission. 

Given  that  the  guidance  will  use  a  reference  control  profile,  the  type  of  reference 
control  profile  must  be  determined.  The  bank  angle  could  be  varied  with  time,  but 
as  explained  in  Section  3.2,  a  bank  profile  which  varies  linearly  with  velocity  provides 
near  maximum  crossrange  for  a  given  downrange.  This  is  the  control  profile  which 
was  used  in  [10]  and  chosen  for  this  guidance.  An  added  benefit  of  this  profile  is  that 
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the  vehicle  flies  near  full  lift  up  at  the  end  of  the  trajectory,  resulting  in  desirable 
shallow  final  flight  path  angles. 

Crossrange  control  is  achieved  by  using  the  lateral  component  of  lift.  The  magni¬ 
tude  of  the  lateral  component  of  lift  is  s;  nply  a  by-product  of  the  bank  angle  chosen 
by  the  downrange  control  logic.  However,  it  is  still  possible  to  use  whatever  magni¬ 
tude  of  lateral  lift  is  available  to  meet  the  crossrange  objectives  by  performing  bank 
angle  sign  reversals.  Changing  the  sign  of  the  bank  angle  changes  the  direction  of  the 
lateral  lift  and  allows  the  vehicle  to  control  crossrange  without  placing  demands  on 
the  magnitude  of  the  bank  angle. 

The  guidance  uses  an  azimuth  controller  to  constrain  the  crossrange.  The  con¬ 
troller  determines  the  need  for  a  roll  reversal  based  on  the  current  vehicle  azimuth 
and  azimuth  to  the  target.  A  reversal  is  commanded  when  the  difference  between 
these  values — the  azimuth  error — exceeds  that  which  is  necessary  to  ensure  that  the 
crossrange  constraint  will  not  be  violated  at  the  target.  The  allowable  azimuth  error 
decreases  as  the  vehicle  nears  the  target. 

Since  the  crossrange  controller  uses  lateral  lift  to  satisfy  its  objectives,  it  is  impor¬ 
tant  that  there  always  be  some  lateral  lift  available.  The  linear  bank  profile  causes 
the  vehicle  to  roll  near  full  lift  up  at  the  end  of  the  trajectory  and  eliminates  the  lat¬ 
eral  component  of  lift  at  a  very  critical  time.  In  order  to  ensure  c^ossrange  authority 
throughout  the  trajectory,  the  basic  linear  bank  profile  is  modified  to  provide  a  mini¬ 
mum  bank  angle,  as  depicted  in  Figure  4.1.  The  minimum  bank  angle,  0mm,  ensures 
that  there  is  always  some  lateral  component  of  lift.  This  reduces  the  amount  of  ver- 
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Figure  4.1.  Linear  Bank  Profile  with  Minimum  Bank  Angle 
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tical  lift  which  can  be  used  to  reach  the  downrange  target  by  not  allowing  the  vehicle 
to  fly  full  lift-up,  and  decreases  the  vehicle  maximum  ranging  capability.  However,  a 
minimum  bank  angle  of  15°  results  in  the  loss  of  less  than  4%  of  the  vertical  lift  while 
maintaining  approximately  25%  of  the  available  lift  for  lateral  control.  Therefore, 
downrange  capability  is  only  slightly  affected,  and  the  trade-off  is  advantageous. 

The  final  dynamic  pressure  is  the  last  constraint  for  which  control  is  designed. 
As  was  seen  in  Chapter  3,  final  dynamic  pressure  is  highly  dependent  on  the  vehi¬ 
cle  design.  Decreasing  the  vehicle  ballistic  coefficient  will  reduce  the  “steady-state” 
or  median  value  of  the  final  dynamic  pressure,  and  if  the  vehicle  is  designed  with 
a  sufficiently  small  ballistic  coefficient,  the  final  dynamic  pressure  can  be  controlled 
passively.  However,  because  decreasing  the  ballistic  coefficient  generally  involves  de¬ 
creasing  the  lift-to-drag  ratio  and  vehicle  ranging  capability,  it  is  desirable  to  reduce 
the  ballistic  coefficient  as  little  as  necessary. 

Another  method  of  controlling  the  final  dynamic  pressure,  as  mentioned  earlier 
in  this  chapter,  is  by  varying  the  bank  angle.  However,  ihe  final  dynamic  pressure  is 
relatively  insensitive  to  the  bank  profile.  The  different  bank  profiles  used  to  generate 
Figures  3.9  and  3.10  in  Chapter  3  produced  limited  changes  in  the  final  dynamic 
pressure.  Most  of  the  variation  in  qj  is  due  to  the  phugoid  nature  of  the  dynamic 
pressure  during  entry. 

Because  the  “steady-state”  value  of  final  dynamic  pressure  is  so  highly  vehicle 
dependent  and  relatively  independent  of  the  bank  angle,  significant  changes  cannot 
be  effected  by  varying  the  bank  angle.  However,  smaller  changes  in  qj  could  be 
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accomplished  if  the  magnitude  of  the  oscillations  in  the  phugoid  were  reduced.  If 
the  entry  vehicle  is  designed  to  satisfy  the  dynamic  pressure  constraint  for  most 
trajectories,  the  guidance  could  be  able  to  effect  small  changes  in  the  vertical  lift 
to  satisfy  the  constraint  for  the  remaining  dispersed  cases  without  disturbing  the 
downrange  control. 

In  order  to  do  this,  an  analytic  dynamic  pressure  controller  was  developed  which 
reduces  the  oscillation  by  driving  the  dynamic  pressure  to  a  constant  value.  This 
controller  is  similar  in  design  to  the  altitude  controller  in  [9]  and  the  heat  rate  ~on- 
troller  in  [10].  While  the  downrange  portion  of  the  guidance  determines  the  bank 
angle,  and  thus  the  vertical  component  of  lift,  which  is  necessary  to  reach  the  target 
in  downrange,  the  dynamic  pressure  controller  determines  a  change  in  that  vertical 
component  of  lift  in  order  to  satisfy  the  dynamic  pressure  constraints.  This  small 
change  or  deviation  in  the  vertical  component  of  lift  shows  up  as  a  sinusoid  on  the 
linear  bank  profile  due  to  the  oscillatory  nature  of  the  dynamic  pressure  history. 

The  sign  of  the  commanded  bank  angle,  5^,  is  determined  by  the  azimuth  con¬ 
troller,  and  the  magnitude  of  the  commanded  bank  angle  is  computed  using  inpu. 
from  the  downrange  and  dynamic  pressure  controllers.  The  commanded  bank  angle 
is  computed  as, 

<t>cmd  =  s 6  arccos  [cos  <t>d  +  6  (cos  <p^\  (4.1) 

where  cos  <f>d  's.  the  downrange  lift  fraction — the  fraction  of  the  currently  available  lift 
which  should  be  in-plane  in  order  to  reach  ilie  desired  downrange,  and  8  (cos  <j>$)  is  the 
dynamic  pressure  lift  fraction — the  c.iange  in  the  in-plane  fraction  of  the  currently 
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available  lift  which  is  desired  to  drive  the  dynamic  pressure  to  the  target  value. 
The  dynamic  pressure  controller  determines  S  (cos  <f>q),  while  the  downrange  controller 
computes  cos  4>d- 

4.3  Downrange  Control 

The  input  of  the  downrange  controller  towards  determining  the  commanded  bank 
angle  is  cos  <f>d ,  the  fraction  of  the  currently  available  lift  which  should  be  in-plane  in 
order  to  reach  the  desired  downrange.  As  explained  in  the  previous  section,  downrange 
control  is  provided  by  a  numeric  predictor- corrector  algorithm  assuming  a  linear  bank 
angle  with  velocity  profile.  The  equation  for  the  bank  profile  used  is, 

4>d  =  <f>des  [ - ~ )  +  4>min  (4-2) 

where  vei  and  vj  are  determined  based  on  the  expected  trajectory  and  is  driven 
by  the  vehicle  capabilities,  vj  is  chosen  a  little  smaller  than  the  actual  expected  final 
velocity  to  allow  for  dispersions  in  the  final  velocity. 

The  predictor  and  corrector  functions  work  together  to  determine  the  particular 
bank  profile  to  be  flown.  Since  a  particular  bank  profile  is  defined  by  its  initial  bank 
angle,  <^ea,  a  new  control  profile  can  be  created  by  simply  changing  <f>de3-  The  corrector 
determines  the  change  in  (f>des  that  will  null  the  downrange  error.  This  A<f>jea  is  based 
on  the  currently  predicted  downrange  error  and  the  sensitivity  of  the  downrange 
error  to  changes  in  <j>dea •  The  predictor  flies  out  the  trajectory  using  the  current  bank 
profile  defined  by  <$>de»  in  order  to  determine  the  currently  predicted  downrange  error. 
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The  predictor  also  flies  out  a  trajectory  using  a  bank  profile  perturbed  in  <^ej.  The 
results  of  these  two  predictions  are  used  by  the  corrector  to  determine  the  sensitivity 
of  the  downrange  error  to  changes  in  <f>dea .  The  procedure  that  the  predictor  and 
the  corrector  follow  in  the  process  of  determining  a  new  bank  profile  is  summarized 
below. 

1.  Set  predictor  desired  bank  angle  to  current  desired  bank  angle,  <f>de3 

2.  Set  predictor  initial  state  to  current  navigation  state 

3.  Predict  downrange  error  at  5  km  altitude 

4.  Set  predictor  desired  bank  angle  to  4>de*+  &<j>des 

5.  Repeat  Steps  2  and  3 

6.  Using  results  of  Steps  3  and  5,  compute  the  sensitivity  of  downrange  to  desired 
bank  angle 

7.  Compute  A <t>dta  based  on  the  results  of  Steps  3  and  6 

8.  Determine  the  new  <f>de3  from  the  old  value  and  the  result  of  Step  7 

Once  a  bank  profile  has  been  chosen  by  the  predictor-corrector,  Equation  4.2  and  the 
current  inertial  velocity  are  all  that  is  necessary  to  determine  the  downrange  bank 
angle,  <j>d.  The  downrange  lift  fraction  is  simply  the  cosine  of  this  angle. 

The  downrange  controller  is  executed  once  at  entry  interface  to  set  the  magnitude 
of  the  initial  bank  angle  command.  In  order  to  save  computational  work,  the  controller 
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is  not  executed  again  until  the  aerodynamic  g-load  exceeds  a  trigger  value,  gbegin- 
Although  the  downrange  controller  produces  a  new  value  for  cos  <f>d  each  guidance 
cycle  after  this,  the  predictor-corrector  is  not  executed  this  frequently  due  to  the 
large  control  authority  available  and  the  short  time  between  predictions  over  which 
dispersions  act  to  perturb  the  trajectory.  Early  in  the  entry  when  the  predicted 
trajectories  are  long  and  computationally  intensive,  the  predictor-corrector  execution 
rate  is  scaled  back  to  a  slower  rate,  Raiw.  Later  in  the  entry,  when  more  frequent 
updates  are  desirable,  the  execution  rate  is  increased  to  Rjat.  The  feist  execution 
rate  begins  once  th*»  vehicle  is  below  a  given  altitude,  hcfat,  and  continues  until  the 
downrange  guidance  is  frozen  at  an  altitude  designated  hcjTZ.  The  value  for  each  of 
these  rates  and  altitudes  is  dependent  upon  the  vehicle  and  the  trajectory,  and  is 
parametrically  determined. 

The  operation  of  the  predictor  and  corrector  functions  is  the  subject  of  the  next 
two  sections. 

4.3.  j  Corrector 

The  corrector  determines  the  sensitivity  of  downrange  10  changes  in  the  desired  bank 
angle,  and  based  on  this  sensitivity  and  the  currently  predicted  downrange  error 
determines  the  proper  4>dea.  The  equations  which  are  used  to  implement  this  strategy 
are  derived  in  this  section. 

Expressing  the  downrange  error  in  a  Taylor  series  expansion  of  the  desired  bank 
angle  and  neglecting  second  order  and  higher  terms  yields, 
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A  DRt  =  -t— ^A  <j>dea  (4.3) 

09de» 

In  order  to  null  the  downrange  error,  the  desired  change  in  downrange  error  must  be 
equal  to  the  negative  of  the  predicted  downrange  error, 


ADRe  =  -DRe  (4.4) 

The  above  equations  combine  to  yield  the  incremental  change  in  the  desired  bank 
angle  which  will  null  the  downrange  error: 


A  4>dea  =  —DR' 

The  partial  derivative  in  the  previous  equation  can  be  approximated  by  the  following 
finite  difference  equation, 

dDR'  DR'  —  <f>i icsp'rt)  ~  DR'  (<j>dea  =  descur ) 

dfa  CS  fidespcrt  fid  cscur 

where  DRe{(j>d',  =  <^eiCTir)  is  the  downrange  error  predicted  using  the  current  bank 
profile,  and  DRe  =  (f>de,P'Tt)  is  the  downrange  error  predicted  using  a  perturbed 
bank  angle  profile. 

The  current  and  perturbed  bank  profiles  are  defined  by, 

fidescur  =  $dea  (4.  I  ) 

and, 

fidespert  =  fides  "b  &<Pdes  (4.8) 


dDR' 


d<f>de 


(4.5) 
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where  is  the  desired  bank  angle  from  the  previous  guidance  execution,  and  8<pdc, 
is  a  small  perturbation  in  desired  bank  angle. 

Based  upon  the  results  of  the  predicted  trajectories  and  Equation  4.5,  the  updated 
downrange  guidance  command  is, 

<t>des  =  <f>deJ  +  A  <j>dea  (4.9) 

Near  the  end  of  the  trajectory,  the  downrange  becomes  very  insensitive  to  changes 
in  4>des  and  the  partial  derivative  in  Equation  4.5  becomes  very  small.  Even  for  small 
downrange  errors,  the  required  incremental  change  in  desired  bank  angle  becomes 
unreasonably  large,  and  this  causes  the  guidance  to  become  unstable.  In  order  to 
deal  with  this  instability,  there  are  three  options. 

At  some  point  in  the  trajectory  prior  to  the  development  of  the  instability,  the 
guidance  can  be  frozen.  Once  the  partial  decreases  below  a  specified  value,  the 
predictor-corrector  is  no  longer  executed,  and  the  control  profile  remains  constant. 
This  avoids  the  instability,  but  the  guidance  loses  its  ability  to  adapt  to  winds — 
which  become  a  major  factor  as  the  vehicle  slows  down — and  other  dispersions  in 
the  environment.  Because  the  instability  generally  occurs  with  significant  time  and 
range  remaining  in  the  trajectory,  the  guidance  must  be  frozen  early — allowing  the 
expected  dispersions  to  cause  large  range  errors.  Since  precision  is  critical,  freezing 
the  guidance  is  not  satisfactory. 

It  is  also  possible  to  move  the  target  aim  point  below  the  required  parachute 
deployment  altitude  so  that  the  instability  is  not  encountered.  Since  control  authority 
is  lost  relatively  early  in  the  trajectory,  this  solution  would  require  an  aim  point 
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significantly  displaced  from  the  desired  terminal  location.  This  is  undesirable  because 
the  guidance  only  controls  to  a  final  state  and  not  to  an  intermediate  state,  so  that  the 
actual  final  conditions  will  differ  greatly  from  the  predicted  termination  conditions. 

The  third  and  chosen  option  is  to  limit  the  magnitude  of  the  incremental  change 
in  desired  bank  angle.  Although  the  partial  is  still  allowed  to  get  small,  the  insta¬ 
bility  that  the  shrinking  partial  causes  is  avoided  by  limiting  the  change  in  (j)de3  each 
guidance  cycle.  This  method  allows  the  guidance  to  remain  responsive  to  dispersions 
throughout  the  trajectory.  However,  it  is  not  desirable  to  limit  the  change  in  desired 
bank  angle  until  the  guidance  has  clearly  established  itself  on  the  desired  trajectory. 
Therefore,  the  limiter  is  not  implemented  until  the  vehicle  has  slowed  below  a  speci¬ 
fied  velocity,  V^/im .  The  magnitude  of  the  limit,  A and  the  appropriate  velocity, 
Va are  specific  to  the  vehicle  and  trajectory  and  are  empirically  determined.  Ad¬ 
ditionally,  in  the  case  of  navigation  updates,  it  may  be  desirable  to  allow  large  changes 
in  the  desired  bank  angle  late  in  the  trajectory.  Although  not  implemented  here,  it 
would  be  possible  to  increase  the  size  of  the  incremental  limit  to  allow  such  changes 
based  on  knowledge  that  an  update  has  occurred.  There  could  be  difficulties  with 
this  due  to  the  instability  of  the  predictor-corrector,  but  they  might  be  handled  with 
appropriate  tests  regarding  the  relative  movement  of  the  target. 

Regardless  of  location  in  the  trajectory,  it  is  possible  that  the  downrange  could 
become  completely  insensitive  to  changes  in  the  bank  angle,  and  the  partial  derivative 
in  Equation  4.5  will  become  zero.  This  corresponds  to  a  complete  loss  of  control 
authority.  If  this  does  occur,  the  incremental  change  in  desired  bank  angle  is  simply 
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set  to  zero  and  the  resulting  desired  bank  angle  is  the  same  as  commanded  in  the 
previous  guidance  cycle. 

4.3.2  Predictor 


The  corrector  requires  predicted  downrange  errors,  as  provided  by  the  predictor,  to 
determine  the  new  desired  bank  angle.  The  predictor  is  a  simplified  3-DOF  trajectory 
simulator.  Using  the  current  bank  angle  profile,  it  numerically  integrates  the  vehicle 
equations  of  motion  from  the  current  time  and  state  until  encountering  termination 
conditions.  The  equations  of  motion  to  be  integrated  are, 


dr 

dt 


=  v 


dv 

Tt  =  °“ 


(4.10) 

(4.11) 


where  the  predictor  initial  conditions  are  the  current  navigation  estimates  of  position 
and  velocity, 

r0  =  r(tcur)  (4.12) 

v0  =  v(tcur)  (4.13) 


The  total  acceleration  on  the  vehicle  is  the  sum  of  the  gravitational  and  aerodynamic 
accelerations, 

=  "b  ®aero  (4-14) 

The  gravitational  acceleration,  which  includes  the  effect  of  the  J2  term,  is, 


a9  = 


2  * 9 


(4.15) 
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where, 


3  r2 

*3  =  ?r  +  2  j^2  ((*  ~  5z2)  &  +  22  (4-16) 

Z  —  lT  •  Ipole  (4.17) 

rr  =  -  (4.i8) 

r 

and, 

rpo/e  =  —  (4.i9) 

V 

The  purpose  of  the  prediction  is  to  determine  downrange  capability.  Since  the 
out-of-plane  component  of  lift  does  not  significantly  contribute  to  the  downrange 
capability,  the  acceleration  due  to  the  out-of-plane  component  of  lift  can  be  neglected 
in  the  prediction.  Thus,  the  aerodynamic  accelerations  can  be  written  in  terms  of  the 
drag  and  vertical  lift  components  of  acceleration  as  such, 


^oero  ^drag  ^drag  T  Qvlift  Ivlift 


where, 


CDqS 

® drag  — 


or, 


depending  on  the  method  of  vehicle  definition,  and 


Qviift 


&drag  COS  tfrp 
V 


9  =  2  Pp 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 
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Vrel  —  Vre!  '  Vrel 

(4.25) 

vrei  =v-{Jm  x  r) 

(4.26) 

Hrag  =  ~lvret 

(4.27) 

-  vret 

tv"‘  ~  |iU 

(4.28) 

=  COS  (f)p  ^  &r)  ^  * vTti 

(4.29) 

and  <f)p  is  the  bank  angle  at  any  instant  in  the  trajectory  based  on  the  bank  angle 
profile, 

/  \ 

V  —  Vj 


<f>p  =  4>de 


+  <f>n 


(4.30) 


vei  -  vf/ 

where  uet,  vj,  and  <f>mtn  are  the  same  values  used  to  define  the  actual  bank  profile, 
and  v  is  the  prediction  estimate  of  inertial  velocity. 

Current  estimates  of  the  vehicle  lift-to-drag  ratio,  (jjj')  >  and  the  atmospheric 
density,  pp,  are  determined  as  explained  in  the  subsection,  “Density  and  L/D  Esti¬ 
mation”. 


Predictor  Integration  Step  Size 

The  predictor  uses  two  different  integration  step  sizes.  A  large  time  step,  A above 
an  altitude  of  h^tsm ,  in  the  predicted  trajectory,  and  a  small  step  size,  At,m/,  below 
this  altitude.  The  large  step  size  reduces  computational  intensity  at  a  small  price 
of  reduced  accuracy.  The  difference  in  the  integration  accuracy  using  the  large  and 
small  time  steps  is  not  significant,  and  the  only  reason  a  small  time  step  was  used 
was  to  obtain  homogeneity  in  the  terminal  altitude. 
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Density  and  L/D  Estimation 

The  accuracy  of  the  predictor  depends  on  the  accuracy  of  the  vehicle  and  environment 
models.  Thus,  it  is  important  that  off-nominal  conditions  in  these  models  are  taken 
into  account.  This  section  explains  how  the  density  and  lift-to-drag  ratio  values  used 
in  the  predictor  are  determined. 

The  predictor  uses  a  table  look-up  implementation  of  a  model  atmosphere  to 
determine  the  expected  density  throughout  the  predicted  trajectory.  The  actual  den¬ 
sity  encountered  during  entry  will  vary  from  this  model,  and  to  compensate  for  these 
density  dispersions,  the  density  used  in  the  predictor  is  modelled  as, 

Pp  —  K p  Pnom  (4.31) 

where  pnom  is  the  density  determined  from  the  nominal  atmosphere  model  as  a  func¬ 
tion  of  altitude  and  Kp  is  a  density  scale  factor  which  is  the  ratio  of  the  measured 
density  to  the  expected  nominal  density, 


(4.32) 


The  actual  density  encountered  is  estimated  based  on  the  measured  drag  acceler¬ 
ation,  the  navigation  relative  velocity,  and  the  nominal  ballistic  coefficient  as, 


-  _  2  |  Oe/j-aj  |  (  TTl 


(cDs)^m  (4'33) 

where  the  measured  drag  is  determined  from  the  inertial  measurement  unit  (IMU) 
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sensed  acceleration,  Simu ,  as, 


l^drajl  —  aIMU 


Vrel 

l^reil 


(4.34) 


Because  there  will  be  short  term  variations  in  the  density,  and  other  error  sources 
which  will  cause  noise  in  the  measured  drag  term,  the  density  scale  factor  must  be 
filtered,  and  the  result  of  Equation  4.32  is  not  actually  used.  The  actual  value  used 
in  Equation  4.31  is  the  output  of  a  first-order  low-pass  filter, 


K,+  =  (1  -KJK,-  +  AT 

The  gain  used  in  this  filter  is  dete’  mined  using, 

AT  =  1  -e~(^) 


Pnom 


(4.35) 


(4.36) 


where  A t  is  the  time  between  successive  acceleration  measurements,  and  rp  is  the  time 
constant  '.ected  to  attenuate  the  high  frequency  noise,  yet  allow  sufficient  response 
to  changes  in  the  density.  Also,  for  the  first  pass  through  the  guidance,  the  value  for 
I\p  must  be  initialized  to  1. 

Since  the  equation  for  the  measured  density  uses  the  nominal  ballistic  coefficient 
and  the  navigation  relative  velocity,  errors  in  these  values  appear  as  density  disper¬ 
sions.  Hence,  the  density  estimator  allows  the  predictor,  and  thus  the  guidance,  to 
be  responsive  to  errors  in  ballistic  coefficient  and  relative  velocity,  as  well  as  density 
dispersions. 

The  lift-to-drag  ratio  will  predominantly  be  in  e:  ror  due  to  the  inability  to  prop¬ 
erly  predict  the  pitch  moment  coefficient.  Errors  in  determining  the  pitch  moment 
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coefficient,  Cmo,  cause  the  vehicle  to  trim  at  an  off-nominal  angle  of  attack  which 
results  in  an  off-nominal  L/D.  The  error  in  L/D  induced  by  Cmo,  as  well  as  the  errors 
in  simply  predicting  Cl  and  Co,  are  the  drivers  behind  the  need  to  estimate  the 
vehicle  L/D  for  use  in  the  prediction. 

The  lift-to-drag  ratio  used  in  the  predictor  is  modeled  as, 


41  =  £ 


DJ, 


o  D, 


1.37 


where  —  is  the  expected  lift-to-drag  ratio  at  the  current  time,  and  Km.  is  the  ratio 

L)  nom  D 

of  sensed  L/D  to  expected  L/D, 

L 

(4.38) 

D  nom 


I<i.=-fLL 

D  L 


The  expected  —  is  computed  in  the  prediction  using  a  table  look-up  with  linear 

U  nom 

interpolation  between  points  to  determine  the  vehicle  aerodynamic  coefficients  based 
on  the  navigation  estimate  of  trim  angle  of  attack.  The  actual  L/D  of  the  vehicle  can 
be  estimated  using  the  measured  drag  and  lift  accelerations  as, 


£ 

D 


\adrag\ 


(4-39) 


with  the  measured  accelerations  determined  from  the  IMU  as  in  Equation  4.34  and 
below, 


I  aHJt, 


=  >/? 


<1  IMU  ■  O-IMU  ~  \ddrog\ 


(4.40) 


As  is  the  case  with  the  density  estimator,  there  will  be  short-term  variations  in 
the  acceleration  measurements.  These  high-frequency  variations  must  not  be  incor- 
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porated  into  the  prediction  through  the  L/D  estimation  because  they  would  cause 
erroneous  predictions.  This  noise  can  be  eliminated  using  a  low-pass  filter  similar  to 
that  developed  for  the  density  estimator,  and  this  result  is  used  in  Equation  4.37, 


L_ 

=  (1  -  I<2)  Kf  4-  K2-jP— 

D  nom 

where  the  gain  is  determined  using, 


K ,  =  1  -e 


(4.41) 


(4.42) 


The  time  constant,  r^,  is  parametrically  determined  to  give  satisfactory  response  to 
variations  in  L/D,  and  the  initial  value  of  K  l.  is  1. 

A  note  worth  mentioning  here  is  that  the  vehicle  model  is  constant  throughout 
the  predicted  trajectory  because  the  simulator  used  a  constant  vehicle  model.  In 
actuality,  the  performance  capabilities  of  an  entry  vehicle  vary  with  flight  regime, 
and  the  prediction  would  probably  require  several  different  models  used  along  dif¬ 
ferent  portions  of  the  trajectory  in  order  to  properly  simulate  vehicle  performance. 
The  constant  vehicle  modeling  was  deemed  sufficient  for  the  demonstration  of  guid¬ 
ance  capabilities,  since  the  objective  of  this  study  is  not  to  develop  a  flight-quality 
algorithm. 


Final  State  Error  Computations 

The  downrange  error  in  the  predictor  is  computed  as  the  difference  between  the 
current  surface  range  to  the  target  and  the  currently  predicted  surface  range  to  be 
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flown.  This  is  graphically  illustrated  in  Figure  4.2.  Since  the  predicted  trajectory  is 


Figure  4.2.  Downrange  Error  Definition 

flown  using  only  the  in-plane  or  vertical  component  of  lift,  the  surface  range  to  the 
predicted  final  position  represents  the  vehicle’s  downrange  capability  only — with  no 
crossrange  component. 

In  order  to  determine  the  downrange  error,  the  current  position  (which  is  also 
the  initial  predictor  position),  the  predicted  final  position,  and  the  target  position 
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all  must  be  known  in  a  planet-fixed  coordinate  system.  Relating  these  vectors  in  the 
inertial  frame  in  which  they  are  otherwise  used  is  not  useful  because  of  the  rotation 
of  Mars. 

First,  the  target  position  vector  must  be  determined  in  planet-fixed  coordinates 
from  the  map  coordinates  and  altitude.  Since  the  map  coordinates  are  given  using 
<f>gd,  the  geodetic  latitude,  the  geocentric  latitude,  <j>gc,  must  be  computed  using, 


< pgc  =  arctan 


'tan  (<(>gd) 


(4.43) 


where, 


-trie)' 


(4.44) 


and  the  planet-flattening  term  which  describes  the  eccentricity  of  the  planet  is, 


-fe) 


(4.45) 


Based  on  the  longitude,  A,  and  geocentric  latitude,  a  unit  vector  in  the  direction 


of  the  target  is  computed  as, 


cos(<j>gc)cos(X) 
itgt=  cos  (<f>gc)  sin  (A) 

sin  {4>gC) 


(4.46) 


Then,  the  target  vector  is  computed  from, 


Ttgt  —  rtgt  ^tgt 


(4.47) 


where  the  magnitude  of  the  target  vector  is  the  sum  of  the  distance  to  the  planet 
surface,  r4,  and  the  target  altitude.  The  distance  to  the  planet  surface  at  the  target 
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location  is  based  on  an  ellipsoidal  planet  model,  and  is  computed  as, 

1  fm 


ra  —  T eq 


Jl  -  /« (2  -  /„)  (l  -  (5 


(4.48) 


Next,  the  initial  and  final  predictor  position  vectors  are  transformed  to  the  planet- 
fixed  reference  frame  as  shown  below. 


1  PREDo 


->  PF 
rPRED( 


=  MfF  ?P 


=  MfFl 


PREDf 


(4.49) 

(4.50) 


where  the  transformation  matrix  between  inertial  and  planet-fixed  coordinates  is  de¬ 
termined  using  fep,  the  time  since  some  epoch  when  the  orientation  of  Mars  is  known, 
and  u>m,  the  rotation  rate  of  Mars,  as, 


MfF  = 


cos(u ;mtep)  -sin  (umtep)  0 
sin  (u >mtep)  cos  (wmiep)  0 
0  0  1 


The  three  position  vectors  of  concern  are  then  unitized, 


*PF 

rPRED0 


PF 


PREDo 


-*  PF 
rPRED0 


(4.51) 


(4.52) 


fPF 

rPREDf 


*PF 

T'V 


-*  PF 
rpREDf 


r  pf 

rPREDf 


r>  PF 

-ill- 

-*  PF 


(4.53) 


(4.54) 


109 


and  the  results  are  used  in  the  following  equations  which  determine  the  predicted 
range  angle  and  the  current  target  range  angle. 

9P  =  arccos  (r"EDo  •  r"EDf )  (4.55) 

0tgt  =  arccos  (r"BDo  •  r")  (4.56) 

The  difference  between  the  predicted  and  target  range  angles  is  the  range  angle  error, 
as  shown  below. 

o„r  =  ep-  etgt  (4.57) 

The  downrange  error  is  then  the  product  of  the  range  angle  error  and  the  planet 
equatorial  radius. 

DRe  =  0eTT  req  (4.58) 

This  equation  assumes  a  spherical  Mars.  Although  this  is  not  the  case,  the  equation 
works  well  because  the  error  angle  is  small.  By  the  definition  of  predicted  downrange 
error,  a  negative  error  means  that  the  vehicle  is  predicted  to  hit  short  of  the  target 
and  a  positive  error  is  past  the  target. 

4.4  Crossrange  Control 

The  final  crossrange  is  controlled  by  limiting  the  deviation  of  the  vehicle  azimuth 
from  the  azimuth  to  the  target.  If  the  vehicle  is  exceeding  the  azimuth  error  limit 
and  the  current  4>cmd  is  causing  the  error  to  increase,  a  bank  reversal  is  called  for. 
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The  magnitude  of  the  maximum  allowable  azimuth  error  is  a  function  of  the  vehicle 
velocity  and  is  defined  by  the  azimuth  error  corridor  in  Figure  4.3.  Initially,  the 


Figure  4.3.  Azimuth  Error  Corridor 


allowable  azimuth  error  is  V’errmax-  This  value  is  chosen  to  keep  the  vehicle  in  the 
general  direction  of  the  target,  but  not  to  so  tightly  limit  the  heading  as  to  cause 
excessive  bank  reversals.  As  the  vehicle  nears  the  target,  it  is  necessary  to  reduce  the 
allowable  error.  At  an  inertial  velocity  of  V^errmo,,  the  error  limit  begins  ramping  down 
to  a  smaller  limit,  0 errmmi  which  is  obtained  at  the  velocity  V^errmiB.  These  values 
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which  define  the  azimuth  error  corridor  are  constants  and  are  empirically  determined. 

The  crossrange  controller  is  executed  once  at  entry  interface  to  ensure  that  the 
vehicle  is  initially  turning  in  the  proper  direction,  and  then  waits  until  the  g-loads 
exceed  the  starting  value,  gbegin ,  to  start  running  at  the  continuous  rate  of  1  Hz. 

The  current  azimuth  is  determined  by  first  converting  the  inertial  relative  velocity 
vector  into  a  North-East-Down  coordinate  frame.  This  is  done  using, 


where, 


M?ed  = 


Vr.,NED  =  V„, 

(4.59) 

.jH* xr'  x(  .-r)  fr“ x  i  r,  1 

l^po/c  X  |^po/c  x  r| 

(4.60) 

In  the  North-East-Down  frame,  the  north  and  east  components  of  velocity  are  then 
used  to  find  the  vehicle’s  compass  heading  from  north, 


fpcur  =  arctan 


(  Vr el*ED\ 
\Vr^ED) 


(4.61) 


The  heading  to  the  target  is  found  by  determining  the  vector  normal  to  the  plane 
which  contains  both  the  current  vehicle  position  vector  and  the  target  vector  in  planet- 
fixed  coordinates, 


-•  -♦  PF  v  -» 

r±  =  r  x  rtgt 


where, 

f  PF  =  MfFr 

From  this  normal  vector,  the  azimuth  to  the  target  is  computed  as, 

iptgt  =  arctan 


(4.62) 


(4.63) 


(4.64) 
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The  current  azimuth  error  is  simply  the  difference  between  the  current  vehicle 
azimuth  and  the  azimuth  to  the  target  as  defined  by, 

0*rr  =  'I’cut  ~  *Ptgt  (4-65) 

If  the  magnitude  of  the  error  is  greater  than  180°,  then  the  error  is  in  the  opposite 
direction,  and  the  error  is  corrected  based  on, 

^Perr  =  V'err  ~  360°  sign  (t perT)  (4.66) 

Based  on  the  definition,  a  positive  azimuth  error  for  an  easterly  trajectory  implies  the 
vehicle  will  fly  south  of  the  target.  In  this  case,  a  bank  reversal  will  be  commanded 
only  if  the  azimuth  error  is  greater  than  the  allowable  limit  and  the  currently  com¬ 
manded  bank  angle,  (pcmdt,  is  positive.  This  second  check,  which  determines  whether 
the  current  bank  angle  is  increasing  or  decreasing  the  azimuth  error,  is  used  to  ensure 
only  one  reversal  is  commanded  for  each  excursion  beyond  the  allowable  azimuth 
limit. 

A  bank  reversal  is  commanded  through  the  azimuth  controller  by  multiplying  the 
current  sign  of  the  bank  angle,  S*,  by  -1.  S'*  is  originally  set  to  be  the  sign  of  the 
initial  bank  angle  and  then  switches  sign  whenever  a  bank  reversal  is  needed.  An 
updated  value  for  S*  is  determined  once  per  second  when  the  guidance  is  cycling. 

4.5  Dynamic  Pressure  Control 

The  dynamic  pressure  controller  uses  an  analytic  control  law  to  compute  deviations 
from  the  linear  bank  profile  which  will  satisfy  the  final  dynamic  pressure  constraints. 
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The  output  of  this  controller  is  6  (cos  (f>g)  which  is  a  desired  change  in  the  in-plane 
fraction  of  the  currently  available  lift. 

The  analytic  control  law  reduces  the  oscillation  in  the  dynamic  pressure  by  driving 
the  dynamic  pressure  to  a  desired  rate,  qtgt ,  and  constant  value,  qigt. 

K-  k~ 

6  (cos  <j>g)  =  -~L  (q  -  qtgt)  +  -^-{q-  qtgt)  (4.67) 

This  control  law  feeds  back  dynamic  pressure  and  the  time  rate  of  change  of 
dynamic  pressure  to  determine  the  change  in  the  lift  fraction  necessary  to  provide  a 
second  order  control  response.  Since  a  constant  dynamic  pressure  is  desired, 

=  0  (4.68) 

The  simplified  analytic  dynamic  pressure  control  law  is, 

6  cos  4>q  =  -£q  +-?-(q-  qtgt)  (4.69) 

The  dynamic  pressure  controller  provides  inputs  to  the  commanded  bank  angle  at 
a  rate  of  1  Hz  once  the  aerodynamic  acceleration  exceeds  the  trigger,  gbegin •  However, 
it  is  desirable  to  use  the  dynamic  pressure  controller  only  when  necessary  to  reduce  qj 
because  the  dynamic  pressure  controller  competes  with  the  downrange  controller  in 
determining  the  magnitude  of  the  bank  angle.  Thus,  the  control  law  is  implemented 
only  when  the  predicted  qj  is  near  or  above  the  desired  limit.  More  specifically,  the 
controller  is  implemented  when  the  predicted  final  dynamic  pressure,  qjp,  exceeds 
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some  fraction  of  the  limit,  or  when 


9/p  >  Kqf  qiim  (4.70) 

where  Kq?  is  a  fraction  which  is  empirically  determined  to  ensure  the  controller  is 
implemented  when  needed  despite  errors  in  the  predicted  final  dynamic  pressure. 

The  predicted  final  dynamic  pressure  is  derived  from  the  downrange  control  pre¬ 
dictor  terminal  conditions  as, 


9/P  =  \pv  vleip  (4-71) 

Additionally,  the  controller  is  not  executed  until  the  vehicle  velocity  has  decreased 
below  a  value,  Vjmp,  which  ensures  that  the  vehicle  has  sufficiently  entered  the  atmo¬ 
sphere.  When  used  at  low  dynamic  pressures  the  controller  tries  to  drive  the  dynamic 
pressure  too  quickly  and  causes  downrange  control  to  be  lost. 

4.5.1  Gain  Determination 

In  order  to  obtain  values  for  the  gains  in  the  controller,  and  to  demonstrate  that  the 
dynamic  pressure  response  can  be  described  as  a  second-order  system,  we  start  by 
writing  the  radial  position  equation  of  motion  for  a  point  mass  entry  vehicle, 


~  m  (cos  7  v)2 

mh  =  L„  -\ - mg 


(4.72) 


where  Lv  is  the  vertical  component  of  the  lift  determined  by, 


Lv  =  Cl<]S  cos  <j>  cos  7 


(4-73) 
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Assuming  a  small  flight  path  angle  yields, 


h  —  (~'L9—  cos  $  +  - - g  (4.74) 

m  r 

Since  the  term  cos  <j>  represents  the  portion  of  available  lift  which  is  in  the  vertical 
direction,  it  can  be  written  as  the  sum  of  the  downrange  lift  fraction  and  the  change 
in  the  lift  fraction  desired  by  the  dynamic  pressure  controller  as, 

cos  <j>  =  cos  <j>d  -f  6  (cos  <t>q)  (4-75) 

Combining  Equations  4.74  and  4.75  yields  the  radial  position  equation  of  motion 
containing  the  lift  fractions  from  both  controllers, 

h  =  [cos  <f>d  +  6  (cos  <f>q)\  +  - - g  (4.76) 

m  r 

Since  the  dynamic  pressure  controller  will  cause  some  perturbation  in  altitude, 

h  =  hd  +  Sh  (4.77) 

and  by  differentiating  twice, 

h  =  hd  +  6h  (4.78) 

The  in-plane  equation  of  motion  becomes, 

hd  +  SK  =  jcos  (j>i  _)_  £  (cog  ^  )]  4.  - - g  (4.79) 

m  r 
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Without  the  dynamic  pressure  controller,  the  in-plane  equation  of  motion  would  sim¬ 
ply  be, 

hd  =  —  [cos  4>d]  +  - - g  (4.80) 

m  r 

Subtracting  the  contribution  of  the  downrange  controller,  Equation  4.80,  from  Equa¬ 
tion  4.79  yields  the  perturbations  which  are  due  to  the  dynamic  pressure  controller 
only. 

Sh-  Qdl[6(cos<t> ?)]  =  0  (4.81) 

m 

Substituting  the  control  law,  Equation  4.69,  into  the  above  equation  yields, 

«  -  [^]  [/^  +  I<q  (?  -  ?,„,)]  =  0  (4.82) 

In  order  to  linearize  this  equation  in  altitude,  the  dynamic  pressure  and  dynamic 
pressure  rate  terms  must  be  expressed  linearly  in  terms  of  altitude.  This  can  be 
done  by  taking  a  Taylor  Series  expansion  of  the  density  term  in  the  simple  dynamic 
pressure  equation. 

?=£  +  fa  h*6h  Vrel 2  (4-83^ 

An  expression  for  the  dynamic  pressure  rate  is  obtained  by  assuming  that  the  time 
rate  of  change  of  relative  velocity  is  small  compared  to  the  change  in  density,  and 
differentiating, 

« =  r-1  ml -Sh  (4-84) 
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Substituting  these  expressions  for  dynamic  pressure  and  dynamic  pressure  rate  into 
Equation  4.82,  yields, 


6h  +  KK-8h  +  KKqSh  =  ~I<q  [\vrel2p(h-)  -  qtgt 


where, 


rs  _  (  ClS\  2  &£_ 

K  ~  2  \  m  )  Vre‘  dh 


h * 


(4.85) 


(4.86) 


An  equation  for  the  derivative  of  the  density  with  respect  to  altitude  can  be  formulated 
if  an  exponential  atmospheric  density  is  assumed  of  the  form, 


P-  p,ie 


~(£) 


(4.87) 


The  partial  derivative  is  then, 


(4.88) 


so  that  the  term  K  is  now  expressed  as, 

K  =  P^  e~(^  (4>89) 

Since  the  expression  for  density  is  linearized  in  the  equations  of  motion,  the  gain  I\ 
must  be  reevaluated  each  time  the  controller  is  executed.  It  is  computationally  more 
efficient  to  gr^up  the  constant  terms  in  K  together  into  KCon,t  and  to  compute  this 
only  once.  KComt  is  then  multiplied  by  the  variable  portion  of  K,  Kvary ,  which  is 
evaluated  each  guidance  cycle.  The  new  values  which  form  the  K  term  are, 


A  —  A  const 


K 


vary 


(4.90) 
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where, 


and, 


K, 


const 


Cl  S  p,i 
2  m  h. 


(4.91) 


hvary  —  VTel 


(4.92) 


Equation  4.85  is  a  second-order  constant  coefficient  differential  equation  in  altitude 
of  ..  form, 

Sh  +  2 (ujh  +  un26h  =  f  ( t )  (4.93) 


The  response  of  this  system  can  be  varied  by  changing  the  damping  ratio  and 
the  natural  frequency.  In  order  to  determine  the  parameters  that  give  the  desired 
transient  response,  the  homogeneous  second-order  differential  equation  in  altitude  is 
considered, 

Sh  +  I<I<-Sh  +  KKqSh  =  0  (4.94) 

The  natural  frequency  and  damping  ratio  of  this  homogeneous  system  are, 


=  jKK~q 

(4.95) 

KKk 

*  c 
& 

II 

V/' 

(4.96) 

and  since  it  is  the  gains  that  are  desired  as  a  function  of  the  input  damping  ratio  and 
natural  frequency,  the  purpose  is  better  served  by  rearranging  as, 

Kq  =  Y  <4-97) 

K--  =  ^  (4.98) 
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A  combination  of  damping  ratio  and  natural  frequency  is  empirically  chosen  to  give 
satisfactory  dynamic  pressure  damping. 
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Chapter  5 

Guidance  Performance 


Different  aspects  of  guidance  performance  are  presented  in  this  chapter.  Section  5.1 
analyzes  the  ability  of  the  guidance  to  steer  the  vehicle  to  all  of  the  targets  that  the 
vehicle  is  capable  of  reaching  under  nominal  open-loop  conditions.  The  next  section 
demonstrates  the  ability  of  the  guidance  to  reach  selected  targets  within  the  vehicle 
footprint  in  the  presence  of  dispersions.  Speed  brake  performance  is  discussed  in 
Section  5.3,  and  the  performance  of  the  dynamic  pressure  controller  is  demonstrated 
in  Section  5.4.  The  last  section  examines  the  behavior  of  the  lift-to-drag  ratio  and 
density  filters  used  in  the  guidance  predictions. 

All  test  entry  trajectories  in  this  chapter,  except  for  the  cases  which  are  dispersed 
in  initial  flight  path  angle,  used  the  same  set  of  state  initial  conditions  summarized  in 
Table  3.1.  Vehicle  M,  described  in  Section  3.4,  was  flown  for  the  tests  at  a  33°  angle 
of  attack  which  provided  an  L/D  of  0.85  and  a  ballistic  coefficient  of  283  kg/sq  m. 
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Table  5.1.  Guidance  I-Loads  —  Baseline  Case  (Nominal  Target) 


Sbegin 

= 

0.07  g’s 

hs 

— 

18288  m 

Vei 

= 

3550  m/s 

Vf 

= 

530  m/s 

Pmin 

= 

15° 

f><f>dea 

= 

3° 

Ralw 

— 

15  sec 

Rjat 

= 

5  sec 

hcfst 

= 

30000  m 

kcjrz 

= 

10000  m 

Aplim 

= 

15° 

= 

2500  m/s 

At  fag 

— 

30  sec 

Atami 

= 

2  sec 

= 

10000  m 

hpqUH 

= 

5000  m 

Kp 

= 

0.0488 

/a 

= 

0.0488 

c 

= 

1.0 

wn 

= 

0.05  rad/s 

K-qs 

= 

0.97 

Vdmp 

= 

0  m/s 

Qtgt 

= 

1000  N/sq  m 

Qlim 

= 

1000  N/sq  m 

^errmai 

= 

6° 

Iperrmin 

= 

3° 

^ernriM 

= 

3300  m/s 

V , 

Verrmin 

= 

2200  m/s 

Predictor  Density  Model  =  COSPAR  Northern  Hemisphere  Mean 


5.1  Nominal  Performance 

The  baseline  performance  of  the  guidance  was  determined  in  a  nominal  atmosphere 
without  the  dynamic  pressure  controller.  The  guidance  constants,  or  I-loads,  con¬ 
tained  in  Table  5.1  were  used  for  these  cases. 

As  seen  in  Figure  5.1,  the  guidance  demonstrated  itself  capable  of  reaching  almost 
all  regions  along  the  edge  of  the  open-loop  footprint  that  was  generated  using  the  bank 
profile  described  in  Chapter  3.  The  short  downrange  part  of  the  footprint  resulted 
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Figure  5.1.  Guidance  Coverage  of  Open-Loop  Footprint 
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in  final  dynamic  pressures  higher  than  allowable,  even  though  the  vehicle  was  able  to 
physically  reach  the  target.  This  is  a  result  of  the  trajectory  simply  being  too  short 
for  the  vehicle  to  dissipate  the  necessary  energy.  The  long  downrange  portion  of 
the  footprint  was  physically  unreachable  because  of  the  introduction  in  the  guidance 
of  <Pmin ■  The  minimum  bank  angle  limits  the  amount  of  vertical  lift  in  order  to 
ensure  some  component  of  lateral  lift  at  all  times.  The  trade-off  is  beneficial  because 
the  maintenance  of  crossrange  authority  is  important  while  the  extra  downrange 
capability  lost  is  negligible. 

The  guidance  resulted  in  the  loss  of  no  crossrange  capability.  As  can  be  seen  in 
Figure  5.1,  the  guidance  was  able  to  achieve  the  same  maximum  crossrange  as  was 
possible  flying  open-loop.  Only  small  portions  of  the  open-loop  footprint  are  lost  to 
the  guidance,  and  these  parts  are  at  the  edge  of  the  downrange  capability  which  is 
not  a  major  concern. 

5.2  Dispersed  Performance 

To  demonstrate  the  capability  of  the  guidance  to  successfully  reach  different  portions 
of  the  closed-loop  footprint  despite  dispersions,  four  points  were  chosen  in  the  foot¬ 
print.  The  four  locations  are  characterized  as  nominal,  short  downrange,  long  down- 
range,  and  maximum  crossrange.  The  location  of  these  targets  within  the  closed-loop 
footprint  is  indicated  in  Figure  5.1.  Runs  were  made  without  the  use  of  the  dynamic 
pressure  controller  to  each  of  these  points  using  the  range  of  dispersions  summarized 
in  Table  5.2.  These  dispersions  consist  of  those  discussed  in  Chapter  2  and  an  initial 
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Table  5.2.  Summary  of  Dispersions 


Aerodynamic 

ct 

±10%  CD 

Coefficients 

Ct 

±10%  CL 

Winds 

wooo 

50  m/s  at  0°  (Northerly  Wind) 

W090 

50  m/s  at  90°  (Easterly  Wind) 

W180 

50  m/s  at  180°  (Southerly  Wind) 

W270 

50  m/s  at  270°  (Westerly  Wind) 

Density 

/ 

Revised  COSPAR  Cool-Low  Model 

pwh 

Revised  COSPAR  Warm-High  Model 

pvl 

Viking  1  Entry  Trajectory  Data 

Pv 2 

Viking  2  Entry  Trajectory  Data 

p~ 

Constant  Density  Bias — 50%  Thin 

p+ 

Constant  Density  Bias — 100%  Thick 

pea 

Early  Density  Shear — 25%  Thin 

Pu 

Late  Density  Shear — 10%  Thin 

Flight  Path 

Angle 

7* 

±0.5°  7et 

flight  path  angle  dispersion.  A  known  flight  path  angle  error  is  the  only  compo¬ 
nent  of  the  state  error  at  entry  interface  which  can  cause  significant  dispersions,  and 
variations  of  ±0.5°  are  considered. 

Using  the  same  I-loads  as  in  the  previous  section,  the  guidance  performed  well 
with  the  exception  of  some  extreme  cases,  and  even  then,  the  results  were  generally 
satisfactory  and  not  unexpected.  The  terminal  conditions  for  the  nominal  and  dis- 
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persed  trajectories  to  each  of  the  four  targets  are  contained  in  Tables  5.3  through  5.6 
on  pages  131-134. 

For  the  nominal  target,  the  guidance  performed  very  well.  None  of  the  continuous 
constraints  were  violated  by  any  of  the  dispersed  trajectories.  The  maximum  g-load 
and  heating  rate  experienced  in  any  of  these  trajectories  are  0.55  g’s  and  13.36  W/sq 
cm/s,  while  the  maximum  integrated  heat  load  was  6643.6  W/sq  cm.  As  seen  in 
Table  5.3,  which  contains  the  nominal  and  dispersed  results  for  the  nominal  target, 
the  only  cases  which  had  difficulties  with  the  terminal  conditions  were  the  constant 
density  bias  trajectories.  The  50%  thin  atmosphere  did  not  provide  enough  drag 
to  slow  the  vehicle  and  caused  the  final  dynamic  pressure  (1012  N/sq  m)  to  be  too 
high.  The  100%  thick  atmosphere  bled  too  much  energy  from  the  trajectory  early 
and  caused  the  vehicle  to  reach  the  target  with  a  flight  path  angle  (-24.7°)  that  was 
too  steep.  Both  of  these  values  barely  exceed  their  respective  limits. 

For  the  short  target,  the  guidance  did  not  perform  as  well,  but  the  overall  perfor¬ 
mance  was  still  very  satisfactory.  For  this  target,  the  nominal  and  dispersed  terminal 
conditions  are  contained  in  Table  5.4.  There  were  no  cases  for  which  the  g-load  or 
heating  rate  exceeded  their  limits,  and  the  maximum  values  for  these  parameters 
were  2.28  g’s  and  17.56  W/sq  cm/s.  The  maximum  total  heat  load  was  4841.7  W/sq 
cm.  The  final  dynamic  pressure  limit  was  exceeded  in  two  cases.  The  shallow  initial 
flight  path  angle  and  the  50%  thin  atmosphere,  respectively,  caused  the  final  dynamic 
pressure  to  be  1202.0  and  1152.2  N/sq  m.  This  is  not  unexpected  because  of  the  short 
length  of  the  trajectory.  The  shallow  initial  flight  path  angle  does  not  allow  the  vehi- 
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cle  to  experience  enough  drag  early  in  the  trajectory,  and  the  thin  atmosphere  does 
not  provide  enough  drag  throughout  the  trajectory  to  sufficiently  reduce  the  velocity. 

For  the  same  reasons  as  in  the  nominal  target  case,  the  terminal  flight  path  angle 
limit  is  exceeded  (-24.6°)  in  the  100%  thick  atmosphere.  The  southerly  (0°)  wind, 
which  is  a  crosswind  for  most  of  the  trajectory,  caused  small  misses  in  both  the 
crossrange  (1.10  km)  and  the  downrange  (1.25  km).  The  crossrange  miss  can  be 
directly  attributed  to  the  wind  causing  the  vehicle  to  veer  off  its  intended  track  late 
in  the  trajectory  when  the  lateral  component  of  lift  is  small.  The  downrange  error 
is  also  caused  by  the  crosswind.  Although  this  may  seem  unnatural,  the  predictor- 
corrector  controls  downrange  using  predictions  which  are  based  on  the  navigation 
estimate  of  the  relative  velocity.  The  navigation  estimate  of  relative  velocity  does  not 
account  for  the  wind,  and  this  causes  a  small  error  in  the  prediction.  The  error  in  the 
prediction  is  not  significant  until  late  in  the  trajectory  when  the  vehicle  slows  down. 
The  slower  velocity  (a  significantly  shorter  vector)  and  a  crosswind  (perpendicular 
vector  which  remains  constant)  can  combine  to  cause  a  large  error  in  the  direction  of 
the  relative  velocity. 

The  only  other  dispersion  which  caused  a  problem  for  the  short  target  is  the  Cool- 
Low  atmosphere.  From  Figure  2.2,  it  can  be  seen  that  this  atmosphere  is  continually 
becoming  thicker  relative  to  the  nominal  atmosphere  as  the  altitude  decreases.  Be¬ 
cause  of  this,  the  downrange  is  controlled  based  on  predictions  which  expect  a  thinner 
atmosphere  than  is  actually  encountered.  The  thicker  atmosphere  increases  the  avail¬ 
able  lift,  and  this  is  why  the  vehicle  flew  slightly  (1.06  km)  past  the  target. 
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Short  targets  generally  pose  a  problem  with  the  final  dynamic  pressure  because  of 
the  reduced  drag  along  the  shorter  trajectories.  The  small  range  errors  at  the  short 
target  can  also  be  attributed,  at  least  in  part,  to  the  short  length  of  the  trajectory. 
Everything  happens  much  more  quickly,  and  the  guidance  has  less  time  to  react  to 
the  dispersions.  Late  in  the  trajectory,  there  is  not  enough  time  with  the  limited 
control  authority  to  completely  correct  the  small  errors. 

With  the  exception  of  one  case,  the  guidance  was  also  very  successful  reaching  the 
long  target.  The  maximum  g-load  and  heating  rate  for  all  of  the  dispersed  cases  to 
this  target  were  0.51  g’s  and  12.21  W/sq  cm/s — well  under  the  limits.  In  addition 
the  maximum  integrated  heat  load  was  7846.0  W/sq  cm.  As  seen  in  Table  5.5,  the 
dispersion  which  caused  a  significant  target  miss  is  the  steep  initial  flight  path  angle. 
The  guidance  commanded  maximum  lift-up  immediately,  yet  the  vehicle  was  still 
unable  to  reach  the  target.  This  is  not  a  failure  of  the  guidance,  but  a  limitation  of 
the  vehicle  lifting  capabilities  because  the  long  target  is  not  inside  what  would  be  the 
dispersed  open-loop  footprint.  Although  not  demonstrated  with  the  same  vehicle,  the 
significant  reduction  in  ranging  capability  caused  by  a  steep  entry  flight  path  angle 
can  be  seen  in  Figure  3.6. 

For  the  long  target,  the  northerly  (180°)  wind  caused  a  small  (1.15  km)  downrange 
error  which  is  again  attributable  to  the  error  in  the  navigation  relative  velocity-an 
error  largely  in  direction  as  opposed  to  magnitude.  The  100%  thick  atmosphere  also 
caused  the  now-expected  error  in  terminal  flight  path  angle  (-24.3°). 

For  the  maximum  crossrange  target,  the  guidance  showed  no  major  irregularities. 
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The  nominal  and  dispersed  results  are  shown  in  Table  5.6.  Once  again,  the  maximum 
g-loads  and  heating  rates,  0.55  g’s  and  13.49  W/sq  cm/s,  were  well  within  acceptable 
values,  and  the  maximum  heat  load  was  only  6550.9  W/sq  cm.  The  50%  thin  atmo¬ 
sphere  caused  a  small  violation  of  the  final  dynamic  pressure  constraint  (1006  N/sq 
m).  And,  just  as  it  did  for  the  other  targets,  the  100%  thick  atmosphere  caused  an 
error  in  the  final  flight  path  angle  (-24.8°). 

In  general,  the  guidance  reacted  properly  to  all  of  the  dispersions  and  performed 
very  well.  Most  of  the  constraint  violations  are  very  small.  Since  the  constraints 
are  not  very  definite  at  this  point  in  the  mission  design,  these  small  errors  are  not 
a  cause  for  concern.  The  only  large  miss  of  the  target  was  shown  to  be  a  vehicle 
limitation  and  not  a  guidance  problem.  The  guidance  demonstrated  the  ability  to 
achieve  precision  entry  using  bank  angle  control  only. 

In  order  to  provide  some  insight  into  the  behavior  of  the  guidance,  the  bank  history 
for  all  four  nominal  cases  to  each  of  the  different  targets  is  plotted  against  velocity 
in  Figures  5.2  through  5.5  on  pages  135-138.  The  effect  of  downrange  on  the  bank 
profile  can  easily  be  seen.  The  minimum  downrange  target  requires  less  vertical  lift, 
so  the  bank  profile  in  Figure  5.3  is  much  steeper  than  the  profile  required  to  reach 
the  nominal  target  in  Figure  5.2.  Along  the  same  lines,  the  profile  necessary  to  reach 
the  maximum  downrange  target  in  Figure  5.4  is  much  shallower  than  the  nominal 
profile  because  more  vertical  lift  is  required  to  reach  the  long  target.  As  can  be  seen 
in  Figure  5.5,  the  maximum  crossrange  case  flew  a  bank  profile  with  approximately 
the  same  slope  as  the  nominal  case.  This  is  because  the  maximum  crossrange  target 
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is  nearly  as  fax  downrange  as  the  nominal  target.  The  difference  in  the  two  cases  is 
in  the  sign  of  the  bank  angle. 

The  azimuth  error  as  a  function  of  velocity  is  presented  for  the  undispersed  trajec¬ 
tories  to  each  of  the  four  targets  in  Figures  5.6  through  5.9  on  pages  139-142.  In  each 
of  the  case s,  the  azimuth  error  becomes  very  large  as  the  vehicle  approaches  or  passes 
the  target.  Although  the  error  gets  very  large,  the  crossrange  controller  still  behaves 
properly  and  only  commands  a  reversal  when  the  current  bank  angle  is  causing  the 
azimuth  error  to  increase. 

Because  nearly  full  lift-up  is  required  to  reach  the  maximum  downrange  target, 
the  lateral  component  of  lift  is  relatively  small  for  this  trajectory.  Hence,  as  can  be 
seen  by  comparing  Figures  5.6  and  5.8,  the  azimuth  error  limit  is  exceeded  fewer  times 
in  this  case  than  in  the  nominal  target  case  which  has  a  larger  lateral  lift  component. 
The  azimuth  error  for  the  maximum  crossrange  target  trajectory  is  initially  relatively 
big,  as  shown  in  Figure  5.9.  Once  this  error  is  reduced  to  within  the  limits,  the  value 
deadbands  between  the  limits  just  as  in  the  nominal  case. 

Plots  of  selected  parameters  from  the  nominal  target,  no  dispersion  case  are  con¬ 
tained  in  Figures  5.10  through  5.16  on  pages  143-149.  The  selected  parameters — 
altitude,  velocity,  dynamic  pressure,  heating  rate,  total  crossrange,  total  downrange, 
and  flight  path  angle — are  plotted  against  time. 
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Table  5.3.  Nominal  Target:  Results  of  Dispersed  Cases 


Dispersion 

DRe 

(km) 

CRe 

(km) 

<7/ 

(N/sq  m) 

7/ 

(deg) 

Nominal 

0.09 

-0.48 

897.4 

-15.3 

Cd 

0.01 

0.45 

860.8 

-16.5 

Cd 

0.14 

0.46 

937.9 

-14.1 

ct 

0.09 

0.33 

839.4 

-15.3 

Cl 

-0.15 

0.04 

986.6 

-15.4 

WOOO 

0.07 

-0.38 

916.8 

-15.9 

W090 

0.56 

-0.47 

897.7 

-15.1 

0.10 

-0.42 

884.3 

-14.9 

W270 

-0.03 

-0.43 

901.5 

-15.5 

Pd 

-0.89 

-0.54 

914.2 

-14.1 

pwh 

0.55 

-0.44 

872.4 

-18.1 

p" 

-0.51 

-0.50 

870.1 

-15.6 

P" 

-0.38 

-0.51 

876.8 

-15.4 

P 

0.32 

0.88 

1012.4 

-9.1 

P+ 

-0.08 

0.23 

869.5 

-24.7 

pes 

0.11 

-0.49 

895.7 

-15.4 

PU 

0.03 

0.00 

910.7 

-13.8 

7 «  + 

0.14 

-0.05 

924.3 

-15.5 

7ct 

889.1 

-15.9 

131 


Table  5.4.  Minimum  Downrange  Target:  Results  of  Dispersed  Cases 


Dispersion 

DRe 

(km) 

CRe 

(km) 

BUB 

7/ 

(deg) 

Nominal 

-0.25 

-0.49 

841.7 

-15.6 

d 

0.18 

-0.46 

837.5 

-16.9 

Ch 

-0.13 

-0.53 

862.7 

-15.1 

Cl 

0.27 

0.45 

792.4 

-16.6 

Cl 

0.07 

-0.54 

934.8 

-15.3 

WOOO 

-1.25 

1.10 

840.4 

-16.2 

W090 

0.71 

-0.43 

840.9 

-14.9 

W180 

-0.79 

0.27 

866.8 

-15.9 

W270 

-0.35 

-0.40 

850.0 

-16.1 

Pcl 

-1.06 

-0.54 

851.2 

-13.8 

0.47 

0.15 

870.6 

-19.0 

Pul 

-0.19 

0.44 

822.3 

-15.9 

Pu2 

-0.27 

-0.52 

820.2 

-15.9 

P~ 

-0.58 

0.31 

1202.0 

mm 

P+ 

0.10 

-0.28 

886.3 

-24.6 

Pea 

0.41 

-0.48 

841.5 

-15.5 

PU 

0.00 

0.00 

850.4 

-13.8 

7  e.+ 

0.32 

-0.42 

1152.2 

-13.9 

7ct 

0.18 

-0.27 

922.7 

-16.5 
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Table  5.5.  Maximum  Downrange  Target:  Results  of  Dispersed  Cases 


Dispersion 

DRe 

(km) 

CRe 

(km) 

.  ^ 

(N/sq  m) 

7/ 

(deg) 

Nominal 

0.25 

-0.55 

919.9 

-15.6 

C+D 

876.3 

-17.0 

Cd 

-0.17 

0.07 

963.2 

-14.2 

Cl 

0.03 

0.06 

865.0 

-15.3 

Cl 

0.70 

0.64 

976.3 

-16.1 

WOOO 

0.90 

-0.43 

955.9 

-15.5 

W090 

0.55 

-0.54 

914.5 

-15.8 

W180 

1.15 

-0.17 

912.5 

-14.9 

W270 

-0.07 

-0.48 

922.1 

-15.4 

Pcl 

-0.75 

-0.58 

937.2 

-14.5 

pwh 

0.57 

-0.46 

893.4 

-18.1 

Pvl 

-0.25 

-0.55 

894.7 

-15.7 

Pv 2 

-0.43 

-0.51 

897.3 

-15.7 

■Mi 

-0.16 

0.41 

945.6 

i 

mmm 

0.03 

864.6 

-24.3 

p” 

-0.14 

922.0 

-15.5 

Pu 

0.22 

0.19 

928.2 

-14.1 

7ei  + 

-0.06 

0.10 

904.9 

-15.8 

Tei 

85.76 

4.64 

916.9 

-15.3 
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Table  5.6.  Maximum  Crossrange  Target:  Results  of  Dispersed  Cases 


Dispersion 

DRe 

(km) 

CRe 

(km) 

•  9/ 

(N/sq  m) 

7/ 

(deg) 

Nominal 

0.05 

0.28 

903.6 

-15.3 

ci 

-0.13 

0.41 

861.7 

-16.6 

Cd 

0.24 

0.47 

944.3 

-14.0 

Ct 

0.10 

0.39 

842.8 

-15.2 

Cl 

0.19 

0.36 

970.1 

•15.3 

wooo 

-0.97 

1.02 

859.3 

-15.6 

W090 

0.90 

0.50 

905.1 

-15.1 

W180 

0.11 

-0.18 

929.7 

-15.2 

W270 

0.09 

0.10 

902.5 

-15.5 

Pcl 
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BANK  ANGLE  (DEO) 


Figure  5.3.  Closed-Loop  Bank  History  for  the  Minimum  Downrange  Target  —  No 
Dispersion  Case 
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BANK  ANGLE  (DEG) 


Figure  5.4.  Closed-Loop  Bank  History  for  the  Maximum  Downrange  Target  —  No 
Dispersion  Case 
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BANK  ANGLE  (DEO) 


Figure  5.5.  Closed-Loop  Bank  History  for  the  Maximum  Crossrange  Target  —  No 
Dispersion  Case 
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12 


Figure  5.6.  Closed-Loop  Azimuth  Error  History  for  the  Nominal  Target  —  No 
Dispersion  Case 
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AZIMUTH  ERRC'  (DEO) 


Figure  5.7.  Closed- Loop  Azimuth  Error  History  for  the  Minimum  Downrange  Target 
—  No  Dispersion  Case 
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Figure  5.8.  Closed-Loop  Azimuth  Error  History  for  the  Maximum  Downrange  Tar¬ 
get  —  No  Dispersion  Case 
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AZIMUTH  ERROR  (DEG) 


Figure  5.9.  Closed-Loop  Azimuth  Error  History  for  the  Maximum  Crossrange  Tar¬ 
get  —  No  Dispersion  Case 
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LTITUDE  (H) 


Figure  5.10.  Closed- Loop  Altitude  History  for  the  Nominal  Target  —  No  Dispersion 
Case 
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Figure  5.11.  Closed- Loop  Velocity  History  for  the  Nominal  Target  —  No  Dispersion 
Case 
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DYNAMIC  PRESSURE  CN/SQ  H) 


Figure  5.12.  Closed-Loop  Dynamic  Pressure  History  for  the  Nominal  Target 
Dispersion  Case 
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Figure  5.13.  Closed- Loop  Heating  Rate  History  for  the  Nominal  Target  —  No 
Dispersion  Case 
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DISTANCE  (KM) 


Figure  5.14.  Closed-Loop  Crossrange  History  for  the  Nominal  Target  —  No  Dis¬ 
persion  Case 
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FLIGHT  PRTH  ANGLE  (DEG) 


Figure  5.16.  Closed-Loop  Flight  Path  Angle  History  for  the  Nominal  Target  —  No 
Dispersion  Case 
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5.3  Speed  Brake  Performance 

Since  entry  vehicles  with  a  high  ballistic  coefficient  tend  to  have  high  terminal  ve¬ 
locities  which  guidance  can  not  significantly  reduce,  the  addition  of  another  method 
of  control  to  the  vehicle  was  considered.  The  effect  of  a  speed  brake,  modeled  as  an 
increase  in  the  coefficient  of  drag,  on  qj  was  analyzed.  A  speed  brake  deployed  early 
in  the  trajectory  decreases  the  final  dynamic  pressure  the  same  as  simply  decreasing 
the  vehicle  design  ballistic  coefficient.  But,  since  a  speed  brake  can  only  cause  small 
changes  in  the  effective  Cb ,  the  decrease  in  the  dynamic  pressure  is  limited. 

Despite  being  limited,  perhaps  the  speed  brake  could  serve  a  useful  purpose  by 
cleaning  up  small  final  dynamic  pressure  errors  in  the  dispersed  cases.  To  test  this 
hypothesis,  the  speed  brake  was  deployed  at  various  velocities  in  a  closed-loop  entry 
trajectory,  and  the  effect  of  the  speed  brake  on  qj  was  examined.  The  entry  trajectory 
was  to  the  nominal  target,  and  the  I- loads  in  Table  5.1  were  used. 

Because  of  the  phugoid  nature  of  the  dynamic  pressure,  deployment  at  one  velocity 
was  found  to  cause  an  increase  in  qj,  while  deployment  200  m/s  later  in  the  same 
trajectory  would  cause  it  to  decrease.  Thus,  the  sensitivity  of  final  dynamic  pressure 
to  velocity  of  speed  brake  deployment  wets  found  to  be  highly  nonlinear.  Due  to  its 
limited  capabilities  and  nonlinear  effect  on  dynamic  pressure,  the  speed  brake  does 
not  appear  to  be  a  feasible  means  of  decreasing  the  final  dynamic  pressure. 
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5.4  Dynamic  Pressure  Control  Performance 

Although  Vehicle  M  was  flown  at  a  high  angle  of  attack  so  as  to  provide  a  smaller 
ballistic  coefficient,  in  the  dispersed  guidance  tests,  there  were  some  cases  in  which 
qf  still  exceeded  the  limit.  The  ability  of  the  dynamic  pressure  controller  to  reduce 
these  values  was  tested  using  the  same  I-loads  as  in  Table  5.1  with  the  exception  that 
Vdmp  was  no  longer  zero,  but  set  to  3000  m/s.  The  results  can  be  seen  in  Table  5.7. 


Table  5.7.  Dynamic  Pressure  Controller  Results 


Case 

Target  Dispersion 

Original  Value 
q}  DRe 

(N/sqm)  (km) 

Controlled  Value 
qj  DRe 

(N/sq  m)  (km) 

1012  0.32 

964  0.09 

1202  -0.58 

1040  -0.03 

MIN  DR  7<,+ 

1152  0.32 

950  3.52 

■i'i  f  I'flal 

1006  -0.32 

976  0.38 

Of  the  four  dispersed  cases  in  which  qj  was  originally  exceeded,  the  controller  was 
successful  at  reducing  the  final  dynamic  pressure  below  the  limit  for  three  of  them. 
To  demonstrate  that  the  controller  functioned  as  designed,  the  dynamic  pressure 
histories  for  the  damped  and  undamped  trajectories  to  the  nominal  target  in  a  50% 
thin  atmosphere  are  presented  for  comparison  in  Figure  5.17.  It  can  be  seen  that  the 
controller  damped  out  the  oscillations  which  resulted  in  a  lower  final  value. 

In  the  third  case  in  Table  5.7 — the  minimum  downrange,  shallow  entry  flight  path 
angle  trajectory — the  final  dynamic  pressure  was  reduced  within  limits  at  the  cost 
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Figure  5.17.  Dynamic  Pressure  Response:  Damped  and  Undamped  Trajectories 
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of  violating  the  downrange  constraint.  This  demonstrates  the  limitation  of  using 
only  one  variable  to  control  three  constraints.  In  the  second  and  unsuccessful  case — 
the  minimum  downrange,  thin  atmosphere  trajectory — the  controller  w as  unable  to 
bring  qj  within  limits  because  the  original  final  dynamic  pressure  was  so  high.  The 
final  value  was  reduced  by  over  150  N/sq  m,  but  because  final  dynamic  pressure  is 
relatively  insensitive  to  bank  angle,  the  controller  can  only  do  so  much.  Further  tests 
using  different  vehicles  with  higher  ballistic  coefficients  produced  similar  results.  The 
dynamic  pressure  controller  was  able  to  damp  the  ohugoid,  but  the  “steady-state” 
dynamic  pressure  value,  which  is  determined  by  vehicle  design,  was  not  reduced,  and 
the  1000  N/sq  m  limit  was  still  exceeded. 

The  dynamic  pressure  controller  is  limited  because  it  trades  ranging  performance 
for  dynamic  pressure  control,  and  because  it  is  unable  to  make  large  reductions  in 
qj.  The  inability  to  greatly  reduce  qj  means  that  vehicle  design  is  the  key  to  making 
significant  reductions  without  resorting  to  some  other  method  of  control.  Given  a 
vehicle  design  which  has  a  reduced  ballistic  coefficient,  the  ability  of  the  controller  to 
consistently  hit  the  target  while  limiting  qj  is  called  into  question  by  the  necessary 
trade-off.  Although  the  dynamic  pressure  controller  performs  as  designed,  active 
control  of  the  final  dynamic  pressure  using  bank  ang^e  only  is  not  promising. 

5.5  Estimator  Performance 

The  performance  of  the  density  and  L/D  estimators  in  the  prediction  is  important  for 
ensuring  guidance  accuracy  in  the  presence  of  dispersions.  A  guidance  which  has  good 
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estimates  of  its  environment  and  vehicle  performance  is  able  to  adapt  to  unpredicted 
variations.  Good  estimates  increase  the  guidance  margin  and  expand  the  edges  of  the 
achievable  footprint  in  the  presence  of  dispersions. 

The  performance  of  the  two  estimators  is  demonstrated  using  an  entry  case  with 
+  10%  Cl  dispersions  and  a  Viking  1  atmosphere.  The  density  estimator  response 
using  a  time  constant  of  20  seconds  in  the  low-pass  filter  is  shown  in  Figure  5.18.  The 
figure  shows  the  actual  density  ratio,  Kp,  and  the  predictor  estimate  of  the  density 
ratio  as  a  function  of  time.  Estimation  does  not  begin  until  approximately  400  seconds 
into  the  trajectory  when  the  sensed  acceleration  exceeds  0.07  g’s.  At  this  point  the 
estimated  density  scale  factor  starts  driving  to  the  actual  value.  After  reaching  the 
actual  value,  the  density  estimator  follows  all  of  the  general  trends  in  the  density.  The 
short-term  density  variations  are  filtered  out  resulting  in  a  smooth  density  estimate 
which  accurately  reflects  the  density  in  the  atmosphere. 

The  L/D  estimator  response  in  Figure  5.19  is  from  the  same  dispersed  trajectory, 
again  using  a  20  second  time  constant.  Shortly  after  the  estimation  begins  around  400 
seconds,  the  estimate  reaches  the  actual  value.  Since  the  L/D  dispersion  is  constant, 
the  estimate  does  not  vary  and  reflects  the  exact  value. 

Both  estimators  worked  exactly  as  desired  in  the  previous  example;  however,  winds 
have  an  adverse  effect  on  the  density  and  L/D  estimates.  Figures  5.20  and  5.21  display 
the  time  response  of  the  density  and  L/D  filters  for  the  same  target  and  dispersions, 
but  with  the  addition  of  50  m/s  crosswinds  (0°  Wind  Direction).  As  seen  in  Fig¬ 
ure  5.20,  approximately  halfway  through  the  trajectory,  the  accuracy  of  the  density 
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Figure  5.18.  Time  Response  of  the  Density  Filter:  No  Wind  Case 
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Figure  5.19.  Time  Response  of  the  L/D  Filter:  No  Wind  Case 
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estimate  begins  to  decrease.  The  estimate  deviates  from  the  actual  density  and  alter¬ 
nates  between  being  high  and  low.  The  source  of  this  deviation  is  the  wind,  and  the 
alternation  is  due  to  roll  reversals.  Wind  induces  an  error  in  the  navigation  relative 
velocity  because  the  model  used  to  determine  vrei  does  not  contains  winds.  Since 
the  measured  density  is  derived  from  the  nominal  ballistic  coefficient,  accelerometer 
measurements,  and  vrei,  wind  also  causes  an  error  in  the  density  estimates.  The  error 
in  the  density  estimate  switches  sign  after  a  roll  reversal  because  the  contribution 
of  the  wind  to  the  relative  velocity  changes  with  the  vehicle  attitude.  For  example, 
the  density  estimate  in  Figure  5.20  changes  from  high  to  low  around  1450  seconds 
as  the  vehicle  changes  direction  and  reduces  the  headwind.  The  error  in  the  density 
estimate  does  not  become  noticeable  until  later  in  the  trajectory  because  the  relative 
error  in  vrei  increases  as  the  vehicle  slows  down.  Although  the  density  estimate  error 
grows,  the  impact  of  this  error  on  downrange  performance  is  negligible  as  the  error  is 
never  very  large  and  is  frequently  changing  signs  with  each  roll  reversal. 

The  L/D  estimator  response  in  Figure  5.21  is  for  the  wind  dispersed  trajectory, 
again  using  a  20  second  time  constant.  Since  the  measured  L/D  is  determined  using 
the  drag  and  lift  acceleration  computations  which  utilize  ure/,  the  error  in  the  estimate 
can  also  be  traced  to  errors  in  vre/,  and  thus  to  winds.  The  L/D  estimator  performance 
is  similar  to  the  density  estimator  performance.  Each  roll  reversal  causes  a  change 
in  the  sign  of  the  error  in  the  L/D  estimate,  and  the  average  error  in  the  estimate 
throughout  the  trajectory  is  near  zero.  Also,  the  error  in  the  L/D  estimate  grows  as 
the  velocity  decreases.  However,  the  error  in  the  L/D  estimate  is  more  pronounced 
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and  is  exhibited  earlier  than  the  error  in  density  estimate.  This  can  be  traced  to 
the  fact  that  the  measured  density  is  based  on  the  sensed  drag  acceleration  and  vrei, 
whereas  the  measured  L/D  is  based  on  the  sensed  drag  acceleration  and  the  sensed 
lift  acceleration.  Since  the  measured  lift  acceleration  requires  both  the  magnitude 
and  direction  of  vfe[,  its  error  due  to  winds  is  greater  than  that  simply  present  in  the 
magnitude  of  ure/. 

If  wind  speed  and  direction  could  be  measured  using  an  on-board  sensor  and 
incorporated  into  the  navigation-derived  relative  velocity  so  that  vrei  more  accurately 
reflected  reality,  the  response  of  the  two  estimators  with  wind  could  be  similar  to  the 
response  in  the  no  wind  case. 
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Figure  5.20.  Time  Response  of  the  Density  Filter:  Wind  Case 
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Figure  5.21.  Time  Response  of  the  L/D  Filter:  Wind  Case 
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Chapter  6 


Conclusions 

6.1  Summary  of  Results 

The  goal  of  this  study  has  been  to  examine  the  precision  entry  problem  at  Mars. 
Various  factors  in  entry  vehicle  performance  at  Mars  were  first  investigated  in  order  to 
determine  basic  vehicle  design  requirements.  The  second  part  of  the  study  presented 
a  guidance  algorithm  for  achieving  precision  entry  with  a  mid-L/D  vehicle  using  bank 
angle  control  only. 

The  basic  goals  for  vehicle  design  are  high  lift-to-drag  ratio  and  low  ballistic 
coefficient.  A  high  L/D  was  shown  to  increase  ranging  capability  as  well  as  reduce  the 
final  dynamic  pressure.  While  a  high  L/D  is  desired,  the  marginal  mission  flexibil’ty 
gained  by  increasing  L/D  beyond  1.5  is  small.  An  L/D  of  1.4  is  needed  to  achieve 
crossranges  approaching  1000  km,  although  a  lower  value  is  sufficient  to  reach  an 
in-plane  landing  site. 

The  final  dynamic  pressure  was  found  to  be  almost  linearly  dependent  upon  the 
vehicle  ballistic  coefficient.  For  a  vehicle  with  an  L/D  of  1.0,  a  ballistic  coefficient  less 
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than  300  kg/sq  m  is  necessary  to  keep  the  final  dynamic  pressure  below  1000  N/sq  m. 
Because  of  the  thin  Martian  atmosphere,  the  ballistic  coefficient  which  will  ensure 
that  the  final  dynamic  pressure  limit  is  not  exceeded  at  Mars  is  significantly  less 
than  the  necessary  ballistic  coefficient  for  Earth  entry.  Since  decreasing  the  ballistic 
coefficient  generally  involves  decreasing  the  L/D,  a  trade-off  must  be  made  between 
range  and  final  dynamic  pressure. 

The  same  trade-off  must  be  made  once  a  vehicle  has  been  chosen.  Given  a  vehicle, 
it  is  possible  to  significantly  alter  performance  by  changing  the  trim  angle  of  attack, 
as  this  varies  the  effective  lift-to-drag  ratio  and  ballistic  coefficient.  It  was  shown  that 
a  vehicle  should  not  necessarily  fly  at  the  trim  angle  of  attack  providing  maximum 
L/D  if  final  dynamic  pressure  is  a  concern.  The  ranging  capabilities  must  be  sacrificed 
to  reduce  the  final  dynamic  pressure. 

The  conclusion  drawn  from  the  short  divert  study  was  that  there  is  sufficient 
control  margin  to  divert  an  entry  vehicle  to  reasonable  alternative  sites  if  the  decision 
to  divert  is  made  early  enough.  This  study  was  done  with  a  low  performance  (L/D 
=  0.85)  vehicle,  and  higher  lift-to-drag  ratios  will  increase  the  divert  footprint.  This 
conclusion  must  be  tempered  with  an  understanding  that  the  study  assumed  nominal 
conditions,  and  adverse  density  dispersions  will  reduce  the  divert  capability  of  the 
vehicle. 

Entry  flight  in  the  thin  Martian  atmosphere  with  a  low-to-mid  L/D  vehicle  was 
shown  to  exhibit  a  large  phugoid  response  in  dynamic  pressure  which  makes  continu¬ 
ous  and  final  constraint  control  difficult.  The  entry  guidance  presented  here  steers  the 
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vehicle  based  on  a  reference  bank  profile  which  is  linear  with  velocity.  The  particular 
reference  bank  profile  that  will  allow  the  vehicle  to  reach  the  target  in  downrange  is 
determined  by  a  numeric  predictor-corrector  algorithm.  The  guidance  also  utilizes 
roll  reversals  based  on  the  vehicle  heading  in  order  to  control  the  final  crossrange. 
Additionally,  an  analytic  control  law  is  superimposed  on  the  linear  bank  profile  to 
reduce  the  oscillations  in  the  dynamic  pressure,  and  thus  reduce  the  final  dynamic 
pressure. 

The  guidance  was  shown  to  provide  near-maximum  crossrange  performance  and 
almost  complete  coverage  of  the  possible  footprint.  It  was  tested  in  the  presence  of 
dispersions  and  demonstrated  the  ability  to  handle  them  well.  This  adaptability  is 
achieved  by  using  in-flight  measurements  to  increase  the  accuracy  of  the  predictions. 

Through  the  guidance  testing,  an  entry  vehicle  with  L/D  =  0.85  and  Cb~  282 
kg/sq  m  was  shown  to  have  the  control  authority  necessary  to  reach  the  target  with 
high  precision.  With  this  vehicle,  the  other  major  constraint,  final  dynamic  pressure, 
was  kept  within  limits  in  most  all  cases.  However,  it  must  be  noted  that  if  a  higher 
ballistic  coefficient  vehicle  is  used,  the  guidance  is  able  to  damp  the  oscillations  in 
dynamic  pressure,  but  it  is  unable  to  significantly  reduce  the  final  value.  Therefore, 
the  necessity  for  a  second  method  of  control,  specifically  angle  of  attack  modulation 
to  reduce  the  final  dynamic  pressure,  depends  on  the  vehicle  design.  If  the  vehicle 
design  passively  reduces  the  final  dynamic  pressure  through  the  proper  selection  of 
ballistic  coefficient  and  trim  angle  of  attack,  then  no  other  method  is  necessary.  But, 
if  the  vehicle  design  is  driven  solely  by  factors  other  than  guidance — such  as  weight 
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and  packaging — resulting  in  a  relatively  high  ballistic  coefficient,  bank  angle  control 
alone  is  limited  and  a  second  method  of  control  will  probably  be  required. 

6.2  Future  Research 

Much  more  research  is  necessary  prior  to  the  development  of  a  flight-quality  Mars 
precision  entry  guidance  algorithm.  One  of  the  major  steps  that  must  be  taken 
is  the  definition  of  specific  mission  requirements.  Once  specific  requirements  have 
been  defined,  further  work  should  address  possible  improvements  and  testing  of  the 
proposed  algorithm. 

There  are  several  areas  in  which  guidance  performance  could  be  enhanced.  The 
incorporation  of  wind  measurements  into  the  navigation  estimate  of  relative  velocity 
would  provide  better  estimates  of  the  L/D  and  density,  and  result  in  more  accurate 
predictions.  Better  predictions  would  increase  the  guidance  accuracy,  or  allow  the 
guidance  to  be  run  at  a  slower  rate.  Crossrange  and  dynamic  pressure  control  would 
also  be  strengthened  by  the  wind  measurements  because  relative  velocity  is  used  in 
these  computations,  too.  The  effect  that  such  measurements,  provided  by  an  air  data 
system,  have  on  performance  must  be  analyzed  to  determine  whether  the  expense  of 
adding  the  system  is  worthwhile. 

Additionally,  the  computational  requirements  of  the  algorithm  can  be  decreased. 
The  predictor  step  size  and  the  corrector  execution  rate  should  both  be  optimized  to 
provide  reduced  computational  loads. 

The  analysis  in  this  thesis  could  be  improved  by  using  the  only  recently  available 
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Mars  GRAM  89  atmosphere  [20],  as  it  provides  wind  and  density  shear  models  which 
axe  much  better  than  the  models  used  here.  The  guidance  could  also  be  further 
tested  by  undertaking  a  complete  guidance-navigation  systems  study.  In  this  study, 
the  interaction  of  the  guidance  with  navigation  updates  from  external  aids  such  as 
beacons,  as  well  as  the  sensitivity  of  the  guidance  to  navigation  errors,  would  be 
investigated. 

Prior  to  the  development  of  a  specific  vehicle  or  mission,  work  can  be  done  regard¬ 
ing  the  effectiveness  of  modulating  angle  of  attack  in  reducing  final  dynamic  pressure. 
If  angle  of  attack  can  effectively  control  the  final  dynamic  pressure,  the  vehicle  de¬ 
signers  may  choose  to  provide  a  means  of  varying  angle  of  attack  so  that  the  vehicle 
can  be  built  with  a  high  ballistic  coefficient. 

Additionally,  a  separate,  more  complete,  study  focusing  on  divert  capabilities 
should  be  accomplished.  Using  different  vehicle  configurations,  the  ability  to  divert 
would  be  tested  with  various  targets  in  the  presence  of  dispersions. 
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